1 - _ - rn 

If/  AD-A1Q0  813  AIR  FORCE  INST  OF  TECH  HR I GHT -PATTERSON  AFB 
If  THE  EFFECT  OF  ROTOR  6E0HETRY  ON  THE  HARMONI 
[S  DEC  80  C  STUERKE.  D  H  NORDOUIST 

1  UNCLASSIFIED  AFIT/SE/EE/80D-34 

OH  SCHOO— ETC  F/6  10/2 
PERFORMANCE  OF  STN— ETCCU 

NL 

M 

W 

1 

J 

j  I 

j  I 

ft 

■ 

^■B 

bbi 

^■B 

■■Bl 

w^m 

BM 

IBBB^ 

■ _ 

r 


AFIT/C.E/EE/80D-34 


Accession  P 

NTIS  GRAScI 
DTIC  TAB 
Unannounced 
Justif icati 

or 

□ 

on_ 

By 

_Distr 

Avai 

Dist 

ft 

lbutio 

labili 

Avail 

Spec 

1/ 

ty  Codes 

and/or 

ial 

» 


THE  EFFECT  OF  ROTOR  GEOMETRY 
ON  THE  HARMONIC  PERFORMANCE  OF 
SYNCHRONOUS  GENERATORS, 

THESIS  « 

AFIT/GE/EE/80  D- 34 

Cecil  Stuerke  David  W.  Nordquist 

Captain  USAF  Captain  USAF 


DTIC 


D 


* 


Approved  for  public  release;  distribution  unlimited. 


AFIT/GE/EE/80D- 34 


THE  EFFECTS  OF  ROTOR  GEOMETRY 
ON  THE  HARMONIC  PERFORMANCE  OF 
SYNCHRONOUS  GENERATORS 


THESIS 

Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  Training  Command 
in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science 


by 


Cecil  Stuerke 
Captain 


Approved 


B.S. 

USAF 

Graduate 

for  public 


David  W.  Nordquist, 
Capta  i  n 

Electrical  Engineering 
December  1980 

release;  distribution  unlimited. 


B.S. 

USAF 


Preface 


In  the  Air  Force,  there  is  a  strong  need  for  the  lightest 
weight  electrical  generator  with  the  highest  possible  efficiency. 
This  brought  about  a  thesis  topic.  In  the  topic  it  was  proposed 
that  two  students  investigate  the  effects  of  rotor  geometry 
on  improving  the  harmonic  performance  of  a  synchronous  gen¬ 
erator.  If  the  harmonic  performance  could  be  improved  enough, 
it  was  hoped  that  wave  shaping  circuitry  external  to  the 
generator  could  be  eliminated. 

We  chose  this  thesis  topic  because  of  the  depth  of  our 
early  graduate  courses  in  power  systems  and  our  keen  desire 
to  increase  our  knowledge  in  this  area. 

We  would  like  to  thank  Dr.  Frederick  Brockhurst  for  his 
guidance  throughout  all  stages  of  this  thesis  effort.  We 
would  also  like  to  thank  Mr.  James  Kummeth  for  his  assistance 
on  the  Finite  Element  analysis  computer  work. 

Cecil  Stuerke 
David  W.  Nordquist 
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Abstract 

The  use  of  rotor  geometry  as  a  specification  in  the 
control  of  generator  voltage  and  flux  density  harmonics  is 
studied.  The  possibility  of  generating  a  specific  voltage 
waveform  at  the  load  with  only  a  three-phase  rectifier  between 
the  AC  generator  and  the  load  is  investigated. 

Since  the  specified  load  voltage  is  DC,  the  internally 
generated  phase  voltage  waveform  must  be  one  which  keeps  con¬ 
stant  the  greatest  instantaneous  potential  difference  between 
any  two  phases. 

The  armature  reaction  MMF  is  expressed  as  an  infinite  set 
of  traveling  waves  (with  respect  to  the  stator)  with  an  infinite 
set  of  velocities.  It  is  shown  that  the  rotor  MMF  cannot  can¬ 
cel  armature  reaction  everywhere  because  the  rotor  MMF  has 
only  the  velocity  of  the  rotor  (not  an  infinite  set  of  velo¬ 
cities)  . 

The  MMFs  present  at  all  the  stator  slots  of  an  armature 
coil  group  are  reflected  into  a  single  slot  in  an  attempt  to 
control  the  air  gap  flux  density  harmonics  by  selectively 
positioning  increments  of  rotor  reluctance.  This  method 
fails  because  a  general  solution  for  the  reflected  flux 
densities  (appropriately  synchronized  with  each  other)  could 
not  be  obtained.  This  method  can  be  further  investigated  if 
such  a  general  solution  can  be  found. 
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Specifying  the  air  gap  reluctance  incrementally  as  the 
rotor  moves  under  successive  stator  slots  is  also  tried  as  a 
means  of  harmonic  control.  This  method  proves  unworkable 
due  to  negative  MMF  requirements  and  conflicting  geometry 
specifications . 

Finite  Element  analysis  is  int:oduced  as  a  computer  aid 
to  the  design  of  a  magnetic  structure.  A  specific  computer 
program,  MAGNETIC,  is  introduced.  Another  computer  program 
is  developed  which  generates  input  data  for  MAGNETIC.  This 
input  data  is  for  the  pole  shoe  and  air  gap  regions  of  the 
generator.  The  program  allows  for  a  variable  air  gap  shape. 

Recommendations  are  included  which  may  allow  the 
successful  specification  of  rotor  geometry  as  a  means  of 
harmonic  control,  or  which  may  prove  the  concept  impossible. 
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THE  EFFECTS  OF  ROTOR  GEOMETRY 


ON  THE  HARMONIC  PERFORMANCE  OF 
SYNCHRONOUS  GENERATORS 

I.  Introduction 


Background 

For  more  than  a  decade,  synchronous  AC  generators  have 
been  used  to  supply  power  to  loads  that  require  other  than 
a  pure  sinusoidal  AC  voltage.  The  circuit  configuration  of 
a  typical  example  is  shown  in  Figure  1. 

The  AC  generator,  also  called  an  alternator,  allows 
sliprings  to  replace  troublesome  DC  generator  commutators. 

In  the  case  of  brushless  AC  generators,  sliprings  are  even 
eliminated.  In  high  voltage  applications,  AC  generators 
have  the  additional  advantages  of  a  more  easily  insulated 
stationary  armature  (Ref  6:336). 

In  order  to  supply  a  non-sinusoidal  voltage,  a  conven¬ 
tional  AC  generator  requires  external  wave  shaping  circuitry. 

A  DC  load  voltage  requires  that  the  AC  generator  output  be 
rectified  and  also  requires  a  large  filter.  A  DC  generator 
application  may  use  a  smaller  filter,  but  this  does  not 
outweigh  the  commutator  disadvantage. 

The  relative  simplicity,  reliability,  and  insensitivity 
to  wear,  corrosion,  and  breakage  make  AC  generators  beneficial 
even  in  DC  applications  (Ref  2:76)  and  (Ref  7).  An  AC  gener¬ 
ator  with  a  tailored  voltage  output  adds  to  these  advantages. 
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The  voltage  in  each  AC  generator  coil  varies  in  time, 
but  is  not  automatically  sinusoidal.  The  voltage  wave  shape 
depends  on  the  magnetic  flux  density  field  in  the  generator's 
air  gap.  The  field  is  distorted  due  to  the  effects  on  salient 
pole  geometry,  conductor  slots,  and  armature  reaction.  The 
time  variation  of  the  distorted  voltage  wave  shape  may  be 
resolved  into  the  sum  of  fundamental  and  higher  order  Fourier 
components . 

The  harmonics  of  the  voltage  time  variation  must  be 
controlled.  The  first  few  harmonics  contain  significant 
power  that  the  prime  mover  must  supply.  If  the  power  from 
some  of  these  harmonics  must  be  trimmed  from  the  load,  it  is 
generally  dissipated  as  wasted  heat  in  a  conventional  gener¬ 
ator  and  filter  circuit.  Without  proper  control,  the  voltage 
wave  shape  inside  the  generator  may  have  peaks  that  require 
increased  armature  insulation.  An  undesirably  large  and 
heavy  generator  is  required  to  solve  the  preceding  problems. 

Conventional  generator  design  includes  methods  for  har¬ 
monic  control.  Generator  geometry  is  normally  used  to  sup¬ 
press  the  most  significant  undesired  harmonics.  Specification 
of  the  distribution,  pitch,  and  skew  of  the  armature  coils 
as  well  as  the  rotor  pole  arc  provides  a  direct  and  effective 
means  of  harmonic  control.  If  the  wave  shape  can  be  trimmed 
to  its  fundamental  component,  the  analysis  necessary  to 
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shape  the  waveform  to  the  needs  of  the  load  is  much  easier. 

For  any  waveform,  the  shaping  amounts  to  the  reintroduction 
of  properly  sized  higher  harmonics  at  the  expense  of  the 
fundamental. 

When  the  generator  is  dedicated  to  specific  load  which 
requires  other  than  the  fundamental  frequency  voltage  and 
current,  the  conventional  design  procedure  is  a  disservice. 

If  the  voltage  wave  shape  is  obtained  directly  at  the  gener¬ 
ator  terminals,  it  is  possible  to  meet  the  load  power  require¬ 
ments  with  a  smaller  filter  and  at  the  same  time  reduce 
generator  losses.  Rotor  geometry  control  can  help  provide 
the  correct  voltage  at  the  generator  terminals. 

Problem  and  Scope 

The  purpose  of  this  thesis  effort  is  to  investigate  the 
possibility  of  generating  a  specific  DC  voltage  waveform  at 
the  load  with  only  a  rectifier  between  the  AC  generator  and 
the  load.  Rotor  geometry  modifications  are  to  be  investigated 
as  the  main  means  of  control.  This  control  method  is  applied 
to  the  generation  of  specific  phase  voltages  (by  controlling 
generator  flux  density  harmonics)  that  will  yield  the  desired 
load  voltage  using  only  rectification  to  modify  the  voltage 
wave  shape.  This  effort  reviews  the  theory  by  which  rotor 
geometry  and  the  voltage  wave  shape  are  related.  From  this 
basis,  processes  for  determining  a  satisfactory  rotor  geometry 
are  studied. 
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As sump t ions 

The  attempts  at  rotor  geometry  determination  through 
required  reluctance  do  not  account  for  the  effects  of  finite 
permeability  or  flux  fringing  in  the  generator  magnetic 
structure.  These  effects  arc  considered  negligible,  but  are 
discussed  as  they  are  observed  in  the  finite  element  analysis 
(Section  X).  The  basic  generator  design  and  the  load  para¬ 
meters  used  in  the  development  of  numerical  examples  are 
assumptions  taken  from  the  generator  design  covered  in 
Appendix  I. 

Approach  and  Presentation 

Section  I  is  an  introduction.  Section  II  reviews  the 
theory  through  which  the  generator  magnetic  flux  density 
is  related  to  rotor  geometry.  This  provides  the  basis  for 
harmonic  control. 

In  Section  III,  the  load  voltage  requirements  are  used 
to  select  the  waveform  of  the  phase  voltage  and  current.  The 
specification  of  the  phase  voltage  waveform  determines  the 
required  magnetic  flux  density  versus  time  in  the  generator. 
The  specification  of  the  phase  current  allows  calculation  of 
the  armature  reaction.  Armature  reaction  MMF  (magnetomotive 
force) is  calculated  in  Section  IV. 

Section  V  develops  an  expression  for  the  rotor  MMF. 
Section  VI  combines  the  armature  reaction  and  rotor  MMF  in 
an  expression  for  the  MMF  at  the  slots  of  a  phase  group  ver¬ 
sus  time.  An  example  is  included. 
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Section  VII  discusses  the  limitations  in  the  use  of 
rotor  MMF  to  cancel  certain  harmonics  of  armature  reaction. 

Section  VIII  describes  an  attempt  to  control  the  flux 
density  harmonics  with  the  introduction  of  selectively  posi¬ 
tioned  increments  of  rotor  reluctance.  Reasons  for  the 
inapplicability  of  this  approach  are  discussed. 

Section  IX  describes  the  incremental  specification  of 
rotor  reluctance  as  it  is  rotated  into  position  under  an 
armature  coil.  The  insufficiency  of  this  approach  is  also 
discussed. 

Section  X  looks  at  the  effects  of  finite  permeability 
and  flux  fringing.  Finite  Element  analysis  is  shown  to  be  an 
excellent  tool  in  fine  tuning  generator  design  to  account 
for  these  effects. 

Section  XI  recommends  other  approaches  to  the  solution 
of  the  harmonic  control  problem  for  further  research. 
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1 1 .  Theory 

This  section  reviews  the  theory  through  which  the 
generator  magnetic  flux  density  is  related  to  rotor  geometry. 
This  discussion  provides  a  basis  for  the  control  of  the  flux 
density  harmonics. 

In  electromechanical  energy  conversion,  the  magnetic 
field  is  the  intermediate  stage.  A  rotating  generator  shaft 
causes  magnetic  lines  of  flux  generated  by  the  rotor  coils 
to  cut  the  armature  coils.  A  voltage  is  induced  in  each 
coil  according  to  Faraday's  Law 


V 


coil 


-dsf>  /dt 


(1) 


where 


(/  ..is  the  instantaneous  induced  voltage 
coil  ° 

<J>  is  the  magnetic  flux  in  webers 

t  is  the  time  in  seconds 


The  coil  voltages  combine  as  coil  group  voltage  which 
combine  as  phase  in  voltages.  The  phase  voltages  are  recti¬ 
fied  to  provide  the  load  voltage. 

"At  a  particular  rotating  speed,  the  instan¬ 
taneous  volts  per  conductor  are  proportional  to 
the  air  gap  flux  density  at  the  conductor.  The 
wave  shape  of  the  conductor  voltage  versus  time 
is  therefore  the  same  as  that  of  the  air  gap  flux 
density  vs  distance  around  the  periphery."  (Ref  3: 

Chap  6,6) 

The  phase  voltage  is  shaped  by  both  the  flux  density  and  the 
layout  connection  of  the  armature  coils. 
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The  shape  of  the  flux  density  versus  distance  may  be 
described  as  the  sum  of  harmonic  components.  Each  component 
contributes  to  the  phase  voltage.  The  stator  coil's  pitch, 
distribution,  and  skew  determine  multiplier  factors  which 
modify  the  voltage  (Ref  5:44-55). 

When  the  coil  width  is  less  than  the  width  of  a  rotor 
pole,  its  voltage  is  modified  by  a  pitch  factor. 


where 


sin 


VTT  W  . 


(2) 


K  is  the  vth  harmonic  pitch  factor 
pv 

v  is  the  space  harmonic  number  being  considered 

W  is  the  winding  pitch  in  radians  or  slots 

T  is  the  first  harmonic  rotor  pole  pitch  in 
P  radians  or  slots 

When  the  coils  are  not  concentrically  placed,  the 
voltage  is  modified  by  a  distribution  factor. 


where 


s  in  (v?t/m) 
q  sin(vir/qm) 


(3) 


K,  is  the  vth  harmonic  distribution  factor 
dv 

q  is  the  width  of  a  phase  belt  in  slots 
m  is  the  number  of  internally  generated  phases 
When  the  coils  are  skewed  with  respect  to  the  rotor 
axis,  the  voltage  is  modified  by  a  skew  factor. 
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both  sources  of  MMl:  rotating  with  respect  to  the  stator. 
The  field  produced  by  these  MMT  sources  may  be  controlled 
through  variation  of  the  magnetic  circuit  reluctance. 


Figure  2  One  Pole  of  the  Rotor  and  the  Stator 
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4  =  MMF/R 


(5) 


where 

MMF  is  the  magnetomotive  force  in  ampere  turns 
R  is  the  reluctance  in  henries  ^ 

Reluctance  of  the  magnetic  path  is: 

R  =  y^a  (6) 

where 

p  is  the  permeability  of  the  material  in  henries/ 
meter 

l  is  the  path  length  in  meters 

? 

a  is  the  area  in  meters 

Because  the  permeability  for  air  is  much  lower  than  the 
permeability  of  the  generator  steel,  the  air  gap  usually 
dominates  the  reluctance  and  the  MMF  drop  around  the  magnet¬ 
ic  circuit.  The  air  gap  reluctance  (above  the  rotor)  travels 
with  respect  to  the  stator  as  does  the  MMF  which  forces  the 
flux  through  it. 

If  both  sides  of  Equation  (5)  are  divided  by  the  area  of 
the  flux  tube  (Ref  1:245),  in  which  the  flux  is  found,  it  is 
evident  that  the  flux  density  is  related  to  the  MMF  through 
the  length  of  the  flux  tube. 

B  =  u  MMF/E  (7) 

m 

where 
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2 

B  is  the  magnetic  flux  density  in  webers/meter 
When  the  reluctance  is  due  to  the  air  gap  distance,  this 
becomes 


B  =  a  MMF/5 (x) 


(8) 


where 


<5(x)  is  the  air  gap  distance  in  meters 

x  is  the  distance  along  the  circumference  of  the 
inside  of  the  stator  in  meters  (stator 
coordinates) 

The  chosen  method  for  tailoring  the  magnetic  flux  density 
is  first  to  calculate  the  MMF  variation  at  each  slot,  and  then 
to  adjust  the  reluctance  through  variation  of  the  length  of 
the  air  gap  at  the  slots.  The  magnetic  flux  density  in  turn 
produces  a  phase  voltage  that  becomes  the  desired  voltage  wave 
shape  when  rectified. 
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III.  Phase  Voltage  and  Current  Wave  Shapes 


This  section  presents  the  theory  and  the  selection  of  an 
acceptable  phase  voltage  wave  shape.  The  current  wave  shape 
which  follows  from  the  voltage  is  also  presented.  The  load 
voltage  requirement  is  the  basic  criteria  on  which  the  selec¬ 
tions  are  based. 

A  number  of  different  phase  voltage  wave  shapes  can  pro¬ 
duce  the  desired  load  voltage  wave  shape  after  being  rectified 
Thus,  the  harmonic  content  of  the  generator  magnetic  field 
that  produces  the  voltage  and  current  wave  shapes  is  a  design 
opt  ion . 


Voltage 

The  generator  must  supply  the  required  load  voltage  in¬ 
spite  of  the  generator  transmission  line  impedances.  If,  as 
is  shown  in  Figure  3,  the  load  voltage  is  specified  and  the 
impedances  are  known,  the  required  line  to  line  voltage  wave 
shape  may  be  determined  by  circuit  analysis. 


where 


V 

g 


(Z, 


ZT) 


(9) 


V 

g 


is  the  internally  generated  line  to  line  voltage 
is  the  required  load  voltage 
is  the  load  impedance  in  ohms 
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I'ijjurc  3  Phase  Circuit  Impedance 


is  the  internal  generator  impedance  in  ohms 


'T 


is  the  transmission  line  impedance  in  ohms 


This  analysis  is  to  be  accomplished  harmonic  by  harmonic  from 
the  components  oi'  the  load  voltage.  Remember  that  the  impe¬ 
dances  are  functions  of  frequency.  It  is  the  phase  voltage 
leading  to  the  composite  1/  that  must  be  induced  in  the 
armature  coils  by  the  aggregate  magnetic  field. 

The  shape  of  the  desired  load  voltage  wave  shape  speci¬ 
fies  the  phase  voltage  difference  (between  any  two  external 
generator  phases)  during  the  period  in  which  it  is  applied  to 
the  load.  This  is  because  the  three-phase  rectifier  applies 
the  greatest  instantaneous  potential  difference  between  any 
two  phases  to  the  load.  If  the  load  voltage  is  DC,  the  maximum 
phase  difference  is  required  to  be  constant.  As  the  voltage 
in  one  phase  falls  below  the  necessary  value,  another  phase 
must  rise  to  replace  it. 

Equally  loaded  symmetrical  phases  will  have  identical 
voltage  wave  shapes.  Successive  phases  will  be  evenly  dis¬ 
placed  through  2tt  electrical  radians  of  rotor  rotation. 

Example  phase  voltage  wave  shapes  which  meet  the  1'C  load 
requirements  are  shown  in  figure  4.  Notice  that  in  figure  4c, 
overlapping  voltages  arc  possible.  In  c  cases,  current  is 

supplied  from  more  than  two  phases  at  a  time. 

The  shape  of  the  phase  voltages  as  they  transition  is 
not  specified  by  t  he  load  voltage  requirements.  It  is 
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unnecessary  for  the  phase  voltage  transition  to  encourage 
proper  current  commutation  because  the  current  in  the  load 
is  controlled  by  the  applied  voltage.  The  rate  of  change 
of  the  phase  voltage  may  be  limited  in  order  to  protect  the 
silicon  controlled  rectifiers  and/or  diodes  used  in  the 
three-phase  rectifier. 

Flexibility  in  the  choice  of  the  phase  voltage  wave 
shape  is  useful  because  this  wave  shape  must  be  controlled 
by  a  properly  engineered  magnetic  flux  density  in  the  air 
gap. 

This  flexibility  may  be  exploited  under  rotor  geometry 
control  in  order  to  accomplish  flux  density  control  with  a 
minimum  mean  air  gap.  This  will  decrease  leakage  and  improve 
generator  performance.  Because  the  phase  voltage  determines 
phase  current,  its  shape  also  affects  the  magnetic  field 
shape  resulting  from  the  armature  reaction. 

The  shape  of  the  load  voltage  during  its  application  to 
the  load  and  its  transition  affects  its  harmonic  content. 

The  proper  choice  of  the  phase  voltage  wave  shape  as  it  is 
applied  to  the  three-phase  rectifier  may  improve  generator 
performance  and/or  simplify  calculations.  Once  the  effective 
rotor  length,  the  rotor  tip  velocity,  and  the  phase  voltage 
are  known,  the  required  magnetic  flux  density  is  also  known. 


coil 


Bv£ 


(10) 
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where 


1/  il  is  the  required  coil  group  voltage 

v  is  the  rotor  tip  velocity  in  meters/second 

l  is  the  effective  rotor  length  in  meters 


Voltage  Example 

See  Appendix  I  (Generator  Design)  for  a  reference  on  the 
generator  specifications  used  in  this  example. 

For  the  example  carried  through  this  thesis  effort,  the 
desired  load  voltage  is  constant  (DC).  Although  the  internally 
generated  phase  voltage  is  properly  derived  from  the  desired 
load  voltage  according  to  Equation  (9) ,  this  example  begins 
with  the  specification  of  the  line  to  line  voltage.  The 
value  of  the  line  to  line  voltage  is  chosen  as  796.74  volts. 
This  is  /3"  times  the  RMS  phase  voltage  of  the  sample  sinusoidal 
AC  generator.  It  is  equal  to  the  RMS  voltage  the  generator 
could  supply  to  a  load  if  it  were  wye  connected  and  if  2^.  and 

Z  were  negligible. 

S 

The  generator  phase  voltage  output  comes  from  the  series 
connection  of  four  coil  groups  (two  coil  groups  per  pole  pair 
from  the  generator's  two  pole  pairs).  The  selected  internal 
phase  voltage  wave  shape  is  shown  in  Figure  5.  The  rectifi¬ 
cation  applies  the  maximum  potential  difference  between  any 
two  external  phases  to  the  load.  Considering  that  this 
difference  is  across  two  external  phases  and  considering  the 
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four  coil  groups  in  series  per  external  phase,  the  required 
coil  group  voltage  is: 

V*coii  =  796.74/8  =  99.59  volts  (11) 

where 

* 

V  ■ ,  is  the  peak  value  of  1/  . , 

coil  *  coil 

The  rotor  tip  speed  (v)  is  119.7  meters/second  and  the 
effective  rotor  length  is  .6438  meters.  The  total  flux  den¬ 
sity  which  the  five  conductors  in  each  coil  group  must  cut 
has  the  same  shape  versus  time  as  V  From  Equation  (10), 

the  maximum  flux  density  is: 

* 

V  . 

B  =  — =  1.292  webers/meter  (12) 

max  vc 

where 

? 

B  is  the  maximum  flux  density  in  webers /meter “ 
max  7 

Figure  5  gives  the  voltage  versus  time  (external  to  the 
generator)  for  the  three  phases.  Phase  B  lags  phase  A  inter¬ 
nally  in  the  generator  by  n/3  electrical  radians  (60  electri¬ 
cal  degrees).  Because  phase  B  is  connected  in  reverse,  it 
lags  phase  A  by  an  additional  u  electrical  radians  (180 
electrical  degrees)  for  a  total  lag  of  4 tt / 3  electrical 
radians  (240  electrical  degrees).  The  flux  density  for  each 
phase  has  the  same  shape  as  the  phase  voltage  variation. 

By  Fourier  analysis,  (Ref  14:Chap  2)  the  harmonic  com¬ 
ponents  of  the  desired  total  flux  density  per  phase  are  given  by 
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Figure  5  Time  Variation  of  the  Chosen  Internal  Phase 
Voltage  and  Flux  Density 


B  =  -  l 
*  g  =  l 


Bmax  g  cos  sin 


(13) 


where 

a)  is  the  frequency  in  radians/second 
g  is  the  time  harmonic  number  being  considered 


Current 

The  phase  current  versus  time  is  determined  by  the  shape 
of  the  phase  voltage  waveform  as  well  as  the  circuit  impe¬ 
dances.  Between  current  commutation  periods,  the  phase 
current  is  the  specified  load  voltage  divided  by  the  speci¬ 
fied  load  impedance.  The  value  of  the  phase  current  during 
commutation  may  be  determined  by  the  use  of  harmonic  compon¬ 
ent  analysis  (Ref  12:400).  Because  the  circuit  may  have 
frequency  dependant  impedances,  each  harmonic  component  of 
the  current  may  have  a  different  angle  by  which  it  lags  the 
voltage  of  the  same  harmonic.  The  composite  waveform  accounts 
for  the  current  lag.  Another  approach  is  to  analyze  the  com¬ 
posite  voltage  waveform  directly  (Ref  7). 

Current  Example 

See  Appendix  I  (Generator  Design)  for  a  reference  on  the 
generator  specifications  used  in  this  example. 

Because  it  is  necessary  to  maintain  the  load  power 
at  100  kilowatts,  the  DC  load  current  is  125.51  amperes  (100 
kilowatts  divided  by  460*/T  volts).  There  are  no  parallel 
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current  paths  inside  the  generator.  Therefore,  between 
commutation  periods,  the  phase  current  is  also  125.51  amperes. 

Figure  6  shows  the  circuit  of  Figure  1  during  a  phase  B 
to  a  phase  A  commutation.  Phase  C  is  conducting  throughout 
this  time,  so  no  diodes  are  shown.  It  is  assumed  that  the 
transmission  line  impedance  and  the  generator  resistance  are 
negligible. 

The  phase  voltage  wave  shape  was  given  in  Figure  5.  A 
constant  voltage  (4  •  (/  ^)  is  applied  across  the  inductive 

coils  of  phases  A  and  B.  The  current  in  phase  B  linearly 
increases  from  zero  and  the  current  in  A  linearly  decreases 
toward  zero.  The  generator  phase  inductance  for  this  example 
is  chosen  such  that  the  current  commutation  angle  is  tt/9 
electrical  radians  or  equivalently  20  electrical  degrees. 

The  current  wave  shape  of  the  three  phases  is  given  in 
Figure  7.  Notice  how  they  follow  the  voltages  of  Figure  5. 

The  time  variation  of  the  current  for  the  three  phases  at 
selected  times  is  given  in  Table  1. 

The  harmonic  components  of  the  phase  current  are: 

I  (t)  =  l  (36  Imax/g*2)  [cos  (gir/9)  -  cos  ( g  2  tt  /  9)  ] 
g=l 

+  i  [sin(g2ff/9)  -  sin(g7i/9)]}  sin  g(wt-n/18)  (14) 

where 

I(t)  is  the  phase  current  in  amperes 

I  is  the  maximum  value  of  the  phase  current 

max 
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Table  1  Normalized  Current  in  Six  Internal 
Phases  versus  Time 


TIME 

PHASE 

(wt) 

A 

B 

C 

A’ 

B' 

C’ 

-  5  7T  / 1 5  . 

-1 

1 

0 

1 

-1 

0 

-  4  TT  /  1  5 

-1 

1 

0 

i 

1 

- 1 

0 

-  3  TT  /  1  5 

-1 

1 

0 

1 

-1 

0 

-  2  TT  /  1  5 

-7 
•  / 

1 

-  .3 

.  7 

- 1 

.  3 

-  TV  / 1 5 

- .  1 

1 

-  .  9 

.  1 

- 1 

.  9 

0 

0 

1 

-1 

0 

-1 

1 

ir/15 

0 

1 

-1 

0 

-1 

1 

2  TT  /  1  5 

0 

1 

-1 

0 

-1 

1 

3  7T  /  1  5 

.  3 

.  7 

-1 

-  .3 

-  .  7 

1 

4  TT  /  1  5 

.  9 

.  1 

-1 

-.9 

-  .  1 

1 

5  TT  /  1  5 

1 

0 

-1 

-1 

0 

1 

6  TT  /  1  5 

1 

0 

-1 

- 1 

0 

1 

7  TT  /  1  5 

1 

0 

-1 

-1 

0 

1 

8  TT  /  1  5 

1 

-.3 

-  .  7 

-1 

.  3 

.  7 

9  TT  /  1  5 

1 

-  .  9 

-  .  1 

-1 

.  9 

.  1 

IOtt/15 

1 

-1 

0 

-1 

1 

0 

IV.  Armature  Reaction 


This  section  develops  an  expression  for  the  time  varia¬ 
tion  of  the  MMF  at  each  of  the  stator  slots  due  to  armature 
react  ion. 

As  magnetic  lines  of  flux  inside  a  generator  cut  armature 
coils,  a  voltage  is  induced  at  the  armature  terminals.  When 
a  load  is  connected  to  the  generator,  current  will  flow 
through  the  armature  coils  generating  an  MMF  known  as  armature 
reaction.  Armature  reaction  MMF  opposes  the  MMF  due  to  the 
rotor  field. 

When  the  MMF  from  the  rotor  coils  is  added  to  the  armature 
reaction  MMF,  the  total  MMF  is  found.  The  desired  magnetic 
flux  density  can  be  produced  provided  that  the  correct  air 
gap  reluctance  is  present. 

Single- Phase  Winding 

A  single  coil  from  each  pole  of  a  single-phase  winding 
produces  a  square  MMF  wave  around  the  air  gap  circumference. 

If  the  coils  span  the  full  width  of  a  pole  (T  )  ,  each  coil 
has  a  width  of  -n  electrical  radians.  The  magnitude  of  the 
MMF  wave  will  oscillate  with  the  frequency  of  the  current  in 
the  coils. 

The  MMF  under  the  armature  coils  may  be  expressed  as 
the  product  of  the  time  variation  of  the  coil  current  and 
the  Fourier  components  of  the  space  variation  of  the  square 
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MMF  envelope.  If  the  x  coordinate  (the  distance  along  the 
circumference  of  the  inside  of  the  stator)  is  zero  at  a  posi¬ 
tion  which  corresponds  to  the  center  of  one  of  the  full 
pitch  coils,  the  space  variation  of  the  MMF  per  conductor  is 
given  by: 


MMF (for  one  conductor) =  I 


max  -it  L ,  v 
v  =  1 


1 

£  -sin(vir/2)  cos(vttx/T  )  (15) 

=  i  v  P 


Because  of  the  term  sin(vir/2),  the  spacial  variation  has 
only  odd  harmonics.  Because  of  the  symmetry,  MMF(-x)  =  MMF(x), 
the  coefficient  of  the  sin(virx/T  )  terms  in  Equation  (15)  are 
all  zero  (and  not  shown). 

The  current  in  the  coils  is  the  phase  current  as 
described  in  Section  III. 


I  (t)  =  I 


l  (a  sin  gut  +  b  cos  gut) 


where 

a  are  the  coefficients  of  the  sine  current  harmonics 
g 

b  are  the  coefficients  of  the  cosine  current  harmonics 
g 

Because:  a  sin(wt  +  <J>)  +  b  cos(ut  +  <j>)  =  a'  sin  ut  +  b'  cos  ut 

a  proper  set  of  coefficients  a  and  b  may  be  established 

g  g 

that  account  for  the  current  lag  angle  of  each  harmonic. 

If,  as  is  usually  the  case,  the  negative  portion  of  the  period¬ 
ic  current  wave  is  a  mirror  image  of  the  positive  position, 
only  odd  harmonics  exist  in  the  time  variation  of  the  current. 
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A  generator  phase  may  contain  several  conductors  dis¬ 
tributed  over  one  pole  pitch  in  a  single  coil  group.  Since 
the  current  in  each  conductor  is  identical,  the  resultant 

MMF  mav  be  summed  using  the  d i s t  r  i  but  i on  factor  K,  .  The 

dv 

space  variation  then  becomes: 


MMF (x) 


MMF 

max 


2 

TT 


00 


l 

v  =  1 


—  sin(vn/2) 

V 


K,  COS  (  V  7T  X  /  T  ) 
dv  v  p 


(17) 


where 


MMF  is  the  maximum  value  of  the  MMF 

max 


MMF 


max 


qs  I 


max 


(18) 


where 

s  is  the  number  of  turns  of  a  coil  group  in  a 
single  stator  slot 

Notice  that  the  coordinate  x  is  now  taken  to  be  :ero  in  the 
center  of  a  coil  group. 

The  space  and  time  variation  of  a  single  phase  is  the 
product  of  the  individual  variations. 


MMF (x , t) 


MMF 


max 


<  l  [a  s in  (gut )  +  b  cos(gut)] 

U=1  ^  g 


*  I  l  ^  s in ( vn / 2 )  Kdy  cos(vttx/T  ) 
"  v  =  l  v 


(19) 


where 

•  indicates  a  product 

g  is  the  time  harmonic  number  being  considered 
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Distribution  Example 

See  Appendix  I  (Generator  Design)  for  a  reference  on 
the  generator  specifications  used  in  this  example. 

The  distribution  factor  for  this  example  is: 


s  in  (  vti  /m)  _  s  i n  ( vir / 6 )  .  7  . . 

q  s  in  ( v  tt  /  mq)  S  sin(v7r/30)  ^  “  ' 

Polyphase  Windings 

If  there  is  more  than  one  coil  group  above  each  pole, 
the  phase  of  each  coil  group  must  be  considered  in  calculating 
the  resultant  MMF.  The  current  in  each  coil  group  is  depend¬ 
ent  upon  the  connection  of  the  internal  machine  phases  to  the 
load.  If  there  are  m  internal  phases,  there  will  be  m  possible 
different  time  variations  of  current  in  the  stator.  bach 
current  wave  shape  will  be  equally  distributed  in  time  with 
a  separation  angle  of  2 tt /m  electrical  radians. 

As  the  current  reaches  maximum  in  successive  phases, 
the  maximum  fundamental  MMF  occurs  at  the  center  of  the  suc¬ 
cessive  coil  groups.  Successive  coil  groups  correspond  to 
successive  internal  machine  phases.  The  MMFs  are  therefore 
distributed  Zn/m  electrical  radians  apart. 

The  MMF  due  to  all  of  the  phases  combined  is  the  maximum 
MMF  times  the  product  of  the  time  K(t,r)  and  space  K(x,r) 
variations  for  all  of  the  m  internal  phases. 


K(t,r)  =  l  [a  sin  g  ( w  t  +  r  2  tt  /  m )  +  b  cos  g ( u t  +  r 2 u / m ) ]  (211 
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where 


K(t,r)  is  the  time  MMF  variation 
r  is  the  individual  internal  machine  phase 

oo 

K(x,r)  =  7  [  [  ~  sin(vir/2)  K,  cos  v  +  r2ir/m)  ]  (22) 

v  =  l  p 


where 

K(x,r)  is  the  space  MMF  variation 

m- 1 

MMF(x.t)  =  MMF  l  [ K  (t , r )  •  K(x,r)]  (23) 

lildA 

r  =  0 


Example 

See  Appendix  I  (Generator  Design)  for  a  reference  on  the 
generator  specifications  used  in  this  example. 

The  space  variation  term  in  Equation  (22)  converges  as  it 
is  summed  over  the  integer  v  from  one  to  infinity.  The 
resulting  value  of  the  expression  for  the  six  internal  phases 
(r  =  0  to  r  =  5)  and  for  the  position  of  the  five  slots  in 
phase  C  (Figure  2)  is  given  in  Table  2. 

Table  2  Normalized  Spacial  Component  of 

Armature  Reaction,  K(x,r),  at  Slots 
of  Phase  C  (Figure  2) 


Phase 

B 

C 

B' 

B 

M 

0 

1 

2 

3 

4 

5 

tr/30 

3tt/30 

5  tt  /  3  0 

7  TT  /  3  0 

9  7i  /  3  0 

.5  .5  -.4  -.5  -.5  .4 

.5  .5  -.2  -.5  -.5  .2 

.5  .5  0  -.5  -.5  0 

.5  .5  .2  -.5  -.5  -.2 

.5  .5  .4  -.5  -.5  -.4 
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With  the  current  wave  shape  that  was  established  in 
Section  III,  it  is  possible  to  evaluate  the  armature 
reaction  MMF  at  each  of  the  slots  versus  time.  According  to 
liquation  (23)  and  using  data  from  Table  1  and  Table  2,  the 
products  of  the  time  variation  and  the  space  variation  may 
be  summed  over  the  applicable  phases.  The  number  of  internal 
phases,  m,  is  equal  to  six  in  this  example. 

Notice  that  both  the  time  and  the  space  variation  terms 
for  phase  A'  are  opposite  in  sign  from  the  entries  for  phase 
A.  The  product  of  the  space  and  time  variation  terms  leads 
to  a  magnetic  field  from  A'  which  is  equal  to  the  magnetic 
field  from  A  both  in  magnitude  and  sign.  Likewise  the  phases 
B  and  B'  and  the  phases  C  and  C'  produce  similiar  results. 

From  Figure  2,  a  positive  current  in  phase  B  produces  an  MMF 
under  the  slots  of  phase  C  which  is  opposite  in  direction  to 
the  MMF  produced  by  a  positive  current  in  phase  A. 

The  entries  in  Table  3  (Total  Normalized  Armature  Reaction 
MMF)  were  obtained  by  changing  the  sign  of  the  phase  B  and  B' 
currents  in  Table  1  and  by  multiplying  the  resultant  Table 
1  array  by  the  transpose  of  the  array  in  Table  2.  This 
multiplication  is  the  direct  normalized  interpretation  of 
Equation  (23).  In  order  to  express  the  actual  MMF,  the 
entries  must  be  multiplied  by  both  the  total  number  of  con¬ 
ductors  in  series  per  internal  machine  phase  and  by  the 
current  which  each  conductor  carries.  Put  another  way,  the 
entries  in  Table  3  multiplied  bv  MMF  (Equation  18)  will 
yield  the  actual  MMF. 
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V.  Rotor  MMF 


This  section  develops  an  expression  for  the  MMF  due  to 
the  rotor  coils. 

The  rotor  coils  are  a  significant  source  of  MMF  in  the 
generator.  Therefore,  the  affect  on  the  controlled  magnetic 
circuit  must  be  accounted  for. 

Two  different  effects  from  the  rotor  are  to  be  investi¬ 
gated  in  the  search  for  a  solution  to  the  problem  of  harmonic 
control.  In  Section  VII,  the  approach  is  to  use  the  rotating 
rotor  field  MMF  (rotating  with  respect  to  the  stator)  to 
cancel  certain  rotating  MMF  (again  rotating  with  respect  to 
the  stator)  components  created  by  the  armature  reaction. 

This  procedure  requires  a  properly  designed  MMF  over  the  rotor 
surface . 

In  contrast,  Section  VIII  and  Section  IX  capitalize  on 
the  reluctance  of  the  rotor  for  control.  For  this  case, 
the  assumption  is  made  that  most  of  the  MMF  drop  is  across 
the  air  gap  reluctance  and  that  the  rotor  MMF  may  be  repre¬ 
sented  (first  order)  as  a  square  wave  over  the  rotor  surface. 
With  these  assumptions,  the  development  of  the  rotor  MMF 
expression  and  its  translation  into  stator  coordinates  is 
continued . 

Without  fringing,  the  width  of  each  square  MMF  pulse  is 

equal  to  the  rotor  pole  arc  b  .  Therefore,  the  pitch  factor 

K  ,  may  be  used  to  form  the  MMF  expression, 
pv 
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MMF 

rotor 


— s  Ic  7  —  sin  (vub  / 2t  )  cos  (vttx'/T  )  (24) 

v  r  F  v  P  P  p 


where 


MMF  is  the  rotor  MMF  in  ampere  turns 

rotor  F 


s  is  the  number  of  turns  in  a  rotor  coil 

r 

Ip  is  the  rotor  current  in  amperes 


bp  is  the  rotor  pole  arc  in  radians  or  slots 

x'  is  the  distance  along  the  outside  circum¬ 

ference  of  the  rotor  in  meters  (rotor 
coordinates) 


The  zero  (origin)  of  the  rotor  coordinate  system  is  taken  at 
a  pole  center. 

The  MMF  expression  for  non-constant  MMF  across  the  pole 
face  may  also  be  achieved  by  Fourier  analysis.  In  this  case, 
saturation  and  fringing  are  the  means  of  achieving  the  desired 
effect.  Finite  Element  analysis  (Section  X)  is  an  important 
tool  to  use  when  considering  saturation  and  fringing. 

Each  harmonic  component  of  the  rotor  MMF  travels  with 
respect  to  the  stator  and  has  the  same  velocity  as  the  rotor. 
The  MMF  due  to  the  rotor  can  be  transformed  into  the  stator's 
set  of  polar  coordinates  by  transforming  each  harmonic  compon¬ 
ent  into  a  forward  traveling  wave. 

00 

MMF  =  —  s  Ir  T  —  sin  (vnb  /2t  )  cos  [  (vttx/t  )-vwt] 

rotor  tt  r  F  v  v  p  p  P 

(25) 
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The  MMF  which  drives  a  flux  tube  may  be  moved  around 
that  tube  without  changing  the  flux  density.  Therefore  with 
radial  flux  lines,  the  traveling  square  wave  can  be  projected 
across  the  air  gap  with  the  same  angular  coordinates. 

As  was  true  of  the  armature  MM T  expression,  the  rotor 
MMF  versus  time  can  he  evaluated  at  any  x  by  substitution  of 
that  value.  Rotor  MMF  is  opposed  by  the  armature  reaction 
MMF.  The  MMFs  may  be  algebraically  summed  to  obtain  an 
expression  for  the  total  MMF  along  the  stator  (Section  VI). 

F.xample 

See  Appendix  I  (Generator  Design)  for  a  reference  on  the 
generator  spec i ficat  ions  used  in  this  example. 

The  chosen  pole  arc  is  2  n  /  3  electrical  radians.  There 
are  533  ampere  turns  in  each  rotor  winding.  The  rotor  MMF  is 
a  square  wave  which  travels  with  the  rotor  under  the  armature 
coils.  The  time  variation  of  the  rotor  MMF  under  the  slots 
of  phase  C  (Figure  2)  is  given  in  Table  4. 


Table  4  Rotor  MMF 
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VI.  Total  MMF 


Section  IV  developed  an  expression  for  the  armature 
reaction  MMI;  and  Section  V  developed  an  expression  for  the 
rotor  MMF.  These  MMFs  are  superimposed  in  the  generator's 
air  gap.  F.ach  MMF  produces  a  magnetic  flux  density  with 
space  and  time  variations  that  are  superimposed.  This  section 
develops  an  expression  for  the  total  MMF. 

At  first  glance  it  seems  possible  to  calculate  the  flux 
density  due  to  armature  reaction  and  then  to  adjust  the  rotor 
reluctance  so  the  rotor  MMF  produces  a  cancelling  flux  den¬ 
sity.  This  idea  is  flawed  hecause  a  change  in  the  magnetic 
circuit  also  changes  t he  magnetic  flux  density  driven  by  the 
armature  reaction  MMF.  Therefore,  instead  of  the  superposi¬ 
tion  of  flux  density,  it  is  the  superposition  of  the  MMF 
which  must  be  accomplished.  The  desired  flux  density  may 
then  be  obtained  from  the  total  MMF  by  adjusting  the  air  gap 
reluctance . 

The  total  MMF  in  the  generator  is  the  sum  of  the 
expressions  in  liquation  (23)  and  liquation  (25).  Because  of 
the  direction  defined  for  positive  current  (Figure  2),  the 
armature  reaction  MMF  and  the  rotor  MMF  add. 

As  is  shown  in  Section  VII,  specification  of  rotor  MMF 
is  not  sufficient  to  control  the  generator  MMF  for  all  of  the 
air  gap  for  all  times.  The  armature  reaction  MMF  produces 
both  space  and  time  variations  that  are  not  controlled  by  the 


rotor  MMF.  Substitution  of  the  space  or  the  time  into  the 
expressions  however,  brings  the  terms  of  both  the  armature 
reaction  MMF  and  the  rotor  MMF  to  a  countably  infinite  series. 
The  coordinates  of  the  stator  conductors  are  important  since 
it  is  at  these  points  that  the  time  variation  of  the  flux 
density  creates  the  desired  voltage  wave  shape.  With  this 
substitution,  control  of  the  harmonics  is  approached  in 
Section  VIII  and  Section  IX. 

Example 

See  Appendix  I  (Generator  Design)  for  a  reference  on 
the  generator  specifications  used  in  this  example. 

The  actual  armature  reaction  MMF  is  obtained  by  multi¬ 
plying  the  normalized  armature  reaction  MMF  in  Section  IV 
by  MMF  (Equation  18) .  The  total  MMF  can  be  obtained  by 
summing  the  armature  reaction  MMF  and  the  rotor  MMF.  This  is 
equivalent  to  multiplying  Table  3  by  MMF  and  adding  the 
result  to  Table  4.  The  total  MMF  is  shown  in  Table  5. 
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Table  5  Total  MMF  at  Phase  C  Slots  versus  Time 
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VII.  MMF  Cancellation 


This  section  considers  rotor  MMF  for  the  cancellation  of 
selected  armature  reaction  MMF  harmonic  components.  For  this 
analysis,  it  is  useful  to  express  both  the  armature  reaction 
MMF  and  the  rotor  MMF  as  traveling  waves. 

Armature  Reaction  MMF 

Equation  (23)  is  the  expression  for  the  armature  reaction 
MMF.  Symmetry  about  the  time  origin  causes  the  cosine  terms 
in  K(t,r)  (Equation  21)  to  be  zero.  The  cosine  term  due  to 
asymmetry  is  subsequently  considered.  For  a  given  g  and  v, 
several  constants  from  Equation  (23)  may  be  combined  by  letting: 

A  =  —  qs  I  a  sin  (vtt/2)  K,  /v  (26) 

gv  it  1  max  g  ^  dv 

where 

A  is  a  constant  whose  value  depends  on  g  and  v 

Remembering  that  the  cosine  coefficient  in  Equation  (21)  is 
zero,  the  armature  reaction  MMF  expression  from  Equation  (23) 
can  be  rewritten  as: 

oo  co  m- 1 

MMF  (x , t )  =  l  l  l  A  sin  g(a>t  +  r2ir/m) 
g=l  v=l  r=0  * 

•  cos  v  [  (irx/T  )  +  (rZrr/m)]  (27) 

Using  the  trigonometric  identity 

sin (a)  cos(E)  =  \  sin(a  -  B)  +  4  sin(a  +  P)  (28) 


4  0 


Equation  (27)  may  be  rewritten  as: 

CO  oo 

MM!'  f  X  ,  t )  l  l  iA  • 

0=1  v=l  6 

|sin  [(gwt)  +vTT.x/Tp) 

+  sin  [(gait)  +  (vTTx/Tp)  -  (g+v)^-] 

+  sin  [(goat)  +  (virx/T  )  -  (g  +  v ) ^] 


+  sin 

[  (gait) 

-  (VTTX/Tp)  ) 

2" 

+  s  m 

[  (gait) 

-  (vmx/Tp)  -  (g-v)  — 

+  sin 

[  (gait) 

-  (VTTX/Tp)  -  (g-v)^L 

+  s  in  [  (gwt)  +  (vTix/Tp)  -  (g  +  v)  ^  ]  +  sin  [  (gait)  -  (vrx/Tp)  -  (g  -  v )  (in— ^  —  ] 

(29) 

Equation  (29)  contains  the  positive  direction  (terms  in  the 
right  column  of  the  equation)  and  negative  direction  (terms 
in  the  left  column  of  the  equation)  traveling  waves. 

When  ip  is  a  multiple  of  m, 


2  IT 

sin  0  +  sin  (0  +  ip— )  +  sin  (q 

m 


,  4tt. 
+  ip — ) 
m 


+  sin(G+<p-— — )  =m  sinC 
v  m 

(30a) 


When  tp  is  not  a  multiple  of  m, 


sin  0  +  sin  (0  +  ’P—)  +  sin  (0  +  ip^-)  +  ....  +  s  in  (0+ip-  — — )  =  0 

(30b) 

This  also  holds  true  for  a  cosine  sum  with  equivalent  arguments. 

Because  the  use  of  sines  and  cosines  in  K(t,r)  preserves 
symmetry  about  the  origin,  g  and  v  are  both  odd  integers. 
Therefore,  their  sums  and  differences  must  be  even. 
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With  ip  analogous  to  (g  +  v)  or  (g  -  v)  in  Equations  (30a) 
and  (30b) 


00  CO 


MMh(x,t)  =  l  I  [t  A  sin(gojt+^)  +  j  A  sin  (gwt-^-)  ] 

.  1  _ 1  fa  g  V  1  _ 


g=l  V=1 


VTTX  ■ 

P 
(31) 

where  for  the  first  sine  term  (g  +  v)  is  an  even  multiple  of  m 
and  for  the  second  sine  term  (g-v)  is  an  even  multiple  of  m. 

When  the  cosine  term  in  K(t,r)  is  not  zero,  the  traveling 
waves  may  also  be  found.  When 


B  =  -  qs  I  b  sin  (vtt/2)  K,  /v 
gv  7r  M  max  g  dv 


(32) 


where 


B  is  a  constant  whose  value  depends  on  g  and  v 
gv 


the  cosine  term  in  Equation  (21)  leads  to: 
oo  oo  m- 1 

MMF  (x  ,  t )  =  III  B  cos  g(wt  +  r27r/m)  cos  v  [  (rrx/T  )  +  (r2ir/m)  ] 
g=l  v=l  r=0  gv  p 

(33) 

The  trigonometric  identity 

cos(a)  cos(S)  =  j  cos(a+3)  +  j  cos(a-g)  (34) 


yields : 


CO  oo 


MMF  (x ,  t )  =  l  l  [j  B  cos  (gcut+  +  j  B  cos  (gut  -  -—)] 

o=l  v=l  °  p  ®  p 


VTTX, 

P 
(35) 


where  for  the  first  cosine  term  (g  +v)  is  an  even  multiple  of  m 
and  for  the  second  cosine  term  (g  j  an  even  multiple  of  m. 
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Therefore,  the  complete  expression  for  the  time  variation 
of  the  magnetic  field  (armature  reaction  MMF  only)  referenced 
to  a  fixed  stator  due  to  the  specified  mtB  internal  phase 
armature  current  is: 


MMF ( x , t ) =  l 
g  =  1 


I 

v  =  l 


It  A 


gv 


IT  X 

s  in  (gojt  +  v— )  + 


y  B 


g^ 


,  ^  ix, 

COS  (gujt  +  v— ) 


,  m  ,  -  ,  .  Tix-.  m  n  ,  .  ix-,, 

~  Aev  sin(g“t-v— )  *  y  B p vj  cos  (gut- v— )  ] 

®  p  ®  ‘  p 

(36) 

As  before,  for  the  first  two  terms,  (g  +v)  is  an  even  multiple 
of  m  and  for  the  last  two  terms,  (g  -v)  is  an  even  multiple 
of  m. 


Example 

See  Appendix  I  (Generator  Design)  for  a  reference  on 
the  generator  specifications  used  in  this  example. 

For  m=6,  the  terms  of  Equation  (36)  are  non-cero  when 
(g+v)  or  (g-v)  is  a  multiple  of  6.  Therefore, 

v  =  +  6k  -  g  (37) 


where 

k  is  an  integer  constant  equaling  0,  +1,  +2,  ... 
Table  6  gives  the  space  harmonic  indexes  (values  of  v) 
corresponding  to  the  odd  harmonics  of  the  time  variation  of 
the  armature  reaction  MMF.  A  forward  traveling  wave  is 
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represented  by  sin(guit  -  v— )  .  The  velocity  of  the  traveling 

P 

wave  components  relative  to  that  of  the  rotor  is  equal  to 
the  negative  time  harmonic  index  (-g)  divided  by  the  space 
harmonic  index  (v) . 


Table  6  Space  Harmonic  Indexes  for  Armature 
Reaction  Traveling  Waves 


The  traveling  wave  for  the  rotor  MMF  is  expressed  in 
Equation  (25).  The  harmonic  index  (v)  in  each  component  is 
the  same  for  the  time  and  space  variation. 

The  cancellation  of  two  orthogonal  functions  (Ref  10:188) 

A  sin(at)  +  B  sin(Bt)  =  0  (38) 

requires  that  the  magnitudes  of  a  and  B  are  equal  or  that 
both  A  and  B  are  zero.  This  applies  to  the  cancellation  of 
armature  reaction  MMF  terms  (liquation  56)  and  rotor  MMF  terms 
Equation  (25). 


If  cancellation  is  to  apply  to  all  positions  (x)  and  all 
times  (t),  the  rotor  MMF  can  only  be  used  in  cancellation  of 
armature  reaction  MMF  terms  in  which  g  =  v.  In  other  words, 
rotor  MMF  will  only  cancel  armature  reaction  MMF  that  has  the 
same  relative  velocity  (g/v)  ;  the  velocity  of  the  rotor. 

Specification  of  rotor  MMF  is  not  sufficient  to  control 
the  generator  MMF  for  all  of  the  air  gap  at  all  times.  The 
armature  reaction  MMF  produces  both  space  and  time  variations 
that  are  not  cancelled.  This  limited  cancellation  is  more 
effectively  accomplished  by  the  specification  of  other  para¬ 
meters  such  as  pole  arc  (see  Section  II). 

Another  concern  with  this  approach  is  that  the  reluctance 
function  controls  the  flux  density  according  to  Equation  (7). 
Thus,  even  if  a  certain  MMF  harmonic  is  cancelled,  there  must 
be  protection  from  the  reintroduction  of  that  harmonic  through 
the  variation  of  the  reluctance.  Therefore,  the  use  of  rotor 
MMF  is  not  a  favored  approach  in  the  control  of  generator 
harmonics . 
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VIII.  Reluctance  Harmonic  Control 

This  section  describes  an  attempt  to  control  phase  voltage 
harmonics  by  using  the  reluctance  to  control  the  flux  density 
harmonics . 

Neglecting  saturation  and  considering  the  air  gap  to  be 
dominant,  the  reluctance  along  the  air  gap  is  proportional  to 
the  air  gap  distance.  In  addition,  because  the  air  gap  rotates 
with  the  rotor,  the  reluctance  versus  distance  has  the  same 
shape  as  the  reluctance  versus  time  at  a  fixed  point  on  the 
stator.  Therefore,  although  the  reluctance  requirements  ma> 
be  solved  on  a  time  basis,  it  can  be  directly  related  back  to 
the  air  gap  and  thus  rotor  geometry.  MMF(x,t)  and  MMF(t)  refer 
in  this  section  to  total  MMF. 

Flux  Density  from  MMF 

Just  as  it  is  the  combined  voltage  of  all  the  conductors 
connected  in  series  in  a  coil  group  which  is  specified  as 
phase  voltage,  it  is  the  sum  of  the  flux  densities  cutting 
these  same  conductors  that  must  be  controlled.  The  flux 
density  at  each  stator  slot  is  determined  by  the  MMF  available 
to  force  the  flux  through  the  reluctance  at  that  point. 

The  MMF(x,t)  around  the  stator  is  known.  It  is  the 
superposition  of  the  MMF(x,t)  due  to  the  armature  reaction 
MMF  (Section  IV)  and  the  MMF  due  to  the  rotor  coils 

(Section  V).  The  time  variation  MMF(t)  at  any  location  along 
the  stator  may  be  obtained  by  substituting  that  location's 
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coordinates  into  the  MMF(x,t)  expression.  For  a  given  MMF(t), 
the  value  of  reluctance  determines  the  magnetic  flux  density. 
Therefore,  for  a  given  MMF(x,t),  t he  desired  phase  voltage 
[Section  III)  may  be  had  by  controlling  the  reluctance. 

The  reluctance  at  any  location  along  the  stator  periphery 
is  also  a  function  of  time.  Therefore,  the  total  magnetic 
flux  density  is: 


1  1  MMF.U) 

HTOT=  x£j  Bx(t)  =  Rx  f  t )  •  Area 


(  3 1) ) 
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pe  r i od . 
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B  =  p/K 


(41) 


If  the  value  of  the  proport ional ity  constant  K  is 
changed,  the  magnitude  of  the  flux  density  is  also  changed. 

The  basis  for  this  concept  involves  a  series  of  these  changes 
combined  with  careful  control  of  the  flux  density  harmonics. 

Reluctance  Specification  Concept 

The  magnetic  circuit  reluctance  function  is  periodic.  The 

reluctance  at  any  slot  is  the  same  function  as  the  reluctance 

at  a  previous  slot  except  for  a  time  delay.  If  the  MMF(t) 

expression  for  each  slot  is  shifted  instead,  the  sum  on  the 

right  side  of  Equation  (39)  has  a  common  denominator.  Using 

* 

the  shifted  and  summed  MMF  (t)  reflected  into  the  first  slot, 
the  solution  of  reluctance  versus  time  will  be  simplified. 

A  pulse  width  current  modulation  technique  for  control 
of  line  voltage  harmonics  was  presented  by  Kr i shnamurthy  (Ref 
11:1269).  His  technique  is  reinterpreted  for  the  time  varia¬ 
tion  of  the  magnetic  circuit.  Here  permeance  is  modulated  to 
control  the  magnetic  field. 

The  resulting  flux  density  variation  B(t)  will  consist  of 
a  set  of  symmetrically  shaped  increments  for  each  set  of  har¬ 
monic^)  controlled  (see  Figure  8c).  The  set  of  increment 
centers  for  the  n^ ^  controlled  harmonic  is  given  the  name  Cn 
(see  Figure  8b)  which  corresponds  to  Bn  in  Reference  11. 


48 


(42) 


B  (C  +  At)  =  B  (C  -  At) 
v  n  v v  n  ' 


and 


B  ( 2k-rr  +  C  )  =  -B  (2kir-C  )  =  B  [(2k+l)n-C  ] 
v  n  v  n  n 


where 


-Bv[ (2k+l)**Cn] 
(43) 


k  is  an  integer  constant  equaling  0,  +1,  +2,  ... 
The  product  of  the  permeance  and  the  MMF  will  be  constant  in 
magnitude  over  the  increments,  but  will  change  sign  every  *t 
radians . 

The  coefficient  of  the  n^  harmonic  of  the  flux  density 
is 

8  B  c 

b  =  — ~~  sin  vd  l  sin  vCn  (44) 

n=l 


where 

d  is  the  width  of  an  increment  (maximum  At) 

c  is  the  number  of  harmonics  controlled  by  the 
increments  of  concern 

When  the  C  's  are  chosen  such  that: 
n 


c 

y  sin  vC  =0  (45) 

'• ,  n 

n=l 


bv  can  be  made  equal  to  zero  for  c  different  harmonics.  There 
will  be  4c  pulses  per  2 tt  electrical  radians. 
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Methodology 

The  harmonics  of  the  flux  density  are  controlled  by 
respec i f icat ion  of  the  air  gap  length  under  the  stator  slots 
during  certain  periods  T  .  The  T  periods  correspond  to  the 
non-zero  periods  of  the  modulated  wave  claculated  for  no 
disturbance  of  the  n  previously  controlled  harmonics.  The 
air  gap  length  over  each  period  is  such  that  although  there 
is  variation  in  the  MMF,  the  flux  density  is  constant.  The 
process  is  accomplished  one  harmonic  at  a  time  with  the  funda¬ 
mental  completed  last. 

For  the  first  non- fundamental  harmonic  to  be  controlled, 

,  the  periods  Tq  correspond  to  the  entire  pole  arc. 

Equation  (8)  is  rearranged  to  solve  for  the  length  of  the 
air  gap  over  the  desired  periods. 

u  MMF ( t ) p 

6(t)  =  - 5 - -  (46) 


The  desired  values  for  the  flux  density  are  given  by  rearrange¬ 
ment  of  Equation  (10). 


v  "  vT 


(47) 


For  the  second  harmonic  to  be  controlled,  v ,  the  air  gap 

is  respecificd  over  the  periods  T. .  The  centers,  C  ,  of 

1  1 

those  periods  arc  chosen  in  each  fundamental  frequency  quarter 
cycle  so  that  Equation  (45)  is  satisfied  (c=l).  Therefore, 
there  will  be  no  contribution  to  the  v. 


harmonic  from  the 


respecification.  The  width  of  the  periods  is  extended  to  the 
point  that  they  meet  each  other  or  the  edge  of  the  rotor  pole. 
The  respec i f i ca t ion  (addition  of  an  increment  of  permeance 
over  the  periods  T^)  is  done  so  that  the  v ?  time  harmonics 
of  the  flux  density  are  the  proper  value. 

Equation  (41)  is  used  to  calculate  the  flux  density  com¬ 
ponent  for  v,  above  the  pole  arc  outside  of  T,.  The  remainder 
of  the  necessary  flux  density,  &B^  ,  must  come  from  the  period 
Tj.  The  flux  density  to  be  added  over  the  periods  T,  must  be 
a  constant  according  to  the  modulation  technique. 

Reluctance  is  proportional  to  6 .  If  the  needed  reluctance 
is  6'  and  the  present  reluctance  is  6,  the  reluctance  to  be 
added  is  the  difference  Returning  to  the  distance  domain, 
the  new  air  gap  length  must  be  6'*  5/(5-<$').  The  respecifi¬ 
cation  for  AB  becomes 


^  =  ^-MMF(t)R  (48) 

v 

If  Tn  extends  over  the  boundaries  of  previous  T, Equation 
(48)  will  be  valid  as  long  as  the  proper  5  versus  time  is  used. 
Further  harmonics  are  controlled  by  additional  increments  of 
re luctance . 


cample 


See  Appendix  I  (Generator  Design)  for  a  reference  on  the 
generator  specifications  used  in  this  example. 
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The  third  harmonic  is  chosen  as  the  first  to  be  controlled. 
The  desired  flux  density  is  :ero.  From  liquation  (46),  ;  is 
infinite. 

T^  is  now  determined  for  control  of  the  fifth  harmonic 
(.the  second  to  be  controlled)  without  disturbance  of  the  third 
harmonic.  The  sum 

1 

r 

'L  ,  s  i  n  v  C 
n=  1  n  n 

is  zero  for  C  =  tt/3.  The  centers,  C.  ,  rr / 3  radians  awav  from 
each  multiple  of  tt  radians  are  tt/5,  2  tt  /  5 ,  4m/5,  and  5i/3. 

The  width  of  each  increment  may  be  up  to  tt/5  before  they  touch. 
At  this  full  width,  the  outside  limits  of  the  increments  T ^ 
are  tt/6  to  5  tt  /  6  and  7  tt  /  6  to  1 1  tt  /  6  .  The  merged  increments  are 
2 tt / 3  wide  corresponding  to  the  width  of  the  chosen  2 tt / 3  pole 
arc.  In  other  words,  a  2 tt / 5  pole  arc  constant  air  gap  pro¬ 
duces  no  third  harmonic.  This  is  consistent  with  conventional 
pole  arc  design  (Ref  9:54). 

The  desired  flux  density  for  the  fifth  harmonic  comes  from 
Section  III  and  liquation  (47).  All  of  this  flux  must  come  from 
T^  because  the  rotor  pole  arc  docs  not  extend  past  T^.  From 
Equation  (48)  and  the  fact  that  6  is  infinite,  6',  (the  limit 
of  [6 ' • 6/ (5 -6  '  ) ]  as  5  goes  to  infinity)  may  be  found. 

Next,  T ^  is  determined  for  control  of  the  next  harmonic 
without  disturbance  of  the  third  and  fiftli  harmonics.  Equation 
(45)  is  simultaneously  solved  for  v  =  3  and  v=5  yielding  C,=  .25a-- 
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and  I...  -  .433r.  The  centers  o!  the  15  (t)  pulses  are  ,233m, 

v  2 

.433",  .  5(>(>;i,  1.233m,  ].433h,  1  .  5h(vn  ,  and  1  .?()(_>: 

radians.  The  periods  meet  and  extend  to  the  edge  of  the  pole 
when  their  width  is  .Hum  radians.  The  reluctance  and  flux 
density  waves  built  with  these  periods  are  shown  in  figures 
8a  and  8c.  The  air  gap  length  over  these  periods  is  adjusted 
to  provide  the  necessary  flux  density  not  provided  outside  of 
Id,.  The  flux  density  adjustment  continues  in  this  manner  for 
other  harmonics.  The  reluctance  over  the  final  periods  1 
which  affect  none  of  the  previously  controlled  harmonics  is 
adjusted  for  the  fundamental  component. 

A  graphical  interpretation  helps  clarify  the  process, 
fhe  periods  of  non- infinite  reluctance,  shown  in  figure  8b, 
are  calculated  as  previously  described.  fhe  air  gap  length 
(proportional  to  reluctance)  is  obtained  during  each  period 
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respectively.  The  curved  bound  of  the  shaded  areas  in  Figures 
8e  and  8f  represent  the  integrand  of  Equation  (50)  when  v  =  3  and 
5  respect i ve 1 y . 

The  signed  shaded  areas  in  Figures 8d  through  8g  represent 
the  integrals  from  Equations  (49)  and  (50).  Algebraically,  the 
signed  areas  of  Figures  8d  through  8g  sum  to  zero.  This 
illustrates  that  Equations  (49)  and  (50)  are  zero  for  v  =  3  and 
for  v  =  5.  Thus,  the  harmonic  components  are  zero  for  these 
harmonics  because  of  the  chosen  reluctance  values  and  periods. 
The  illustration  is  similiar  for  other  higher  harmonics. 

How  Many  Harmonics 

The  process  as  described,  approaches  the  desired  flux 
density  wave  shape  one  harmonic  at  a  time.  To  see  the 
advantage  of  this  approach,  consider  specification  of  the 
reluctance  as  a  composite  wave  shape. 

6(1)  =  --B>(t)t'>  (51) 

where 


<5(t)  is  a  time  varying  air  gap 

This  requires  a  specific  finite  reluctance  at  all  times.  When 
the  large  air  gap  between  the  poles  is  under  the  coil  group, 
the  reluctance  is  not  variable.  This  shortcoming  affects  all 
the  harmonics.  Therefore,  the  chosen  method,  because  the 
important  harmonics  can  be  directly  controlled  has  an  advantage. 
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The  number  of  harmonics  which  can  be  controlled  is 
dependent  on  those  chosen  and  the  pole  arc.  Equation  (45) 
yields  Cj  =  .  1 S  5  rr  ,  C,  =  ,344tt,  and  C,  =  ,433tt,  for  control 
of  the  3rd,  5th,  and  7th  harmonics.  If  the  pole  arc  length 
was  less  than  .  633tt,  C,  would  be  at  the  pole  edge  and  control 
of  these  three  harmonics  could  not  be  achieved.  Therefore, 
the  harmonics  and  the  pole  arc  must  be  selected  consistent 
with  the  methodology. 

The  Contradiction 

The  concept  of  harmonic  control  appears  to  be  practical 

once  MMFn  (the  reflected  MMF)  is  obtained  and  the  control 
K 

harmonics  are  selected.  However,  obtaining  MMF^  is  not  an 
easy  problem. 

The  proposed  shift  of  MMF  at  consecutive  slots  1  and  2 
is 

MMF1*(t)  =  MMF  2 ( t - A  t )  (52) 

where 

MMF^  (t)  is  the  shifted  MMF 

The  reluctance  at  one  slot  is  the  same  function  as  the 
reluctance  at  the  previous  slot  except  for  a  delay  of  At. 

In  other  words 

Rj(t)  =  R2(t  -  At)  (53) 

An  expression  for  the  flux  density  at  the  second  slot  is 
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1*2  (t  -  At) 


(54) 


MMFZ (t-At) 


The  shifted  MMF  and  reluctance  are  then  plugged  into 
Equation  (54).  The  contradiction  is  that  the  flux  density 
B^it-At)  reflected  into  the  first  slot  is  not  synchronized 
with  the  flux  density  at  the  first  slot  B^(t).  Such  a  sum 
makes  little  sense. 

A  general  solution  for  a  reflected  B2(t)  synchronized 
with  B^(t)  and  with  the  knowledge  of  MMF?(t),  but  not  R?(t) 
exceeds  the  resources  of  this  thesis  effort.  Therefore, 
another  approach  is  preferred. 
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IX. 


Incremental  Reluctance  Spec i f  i  c a t i on 


This  section  describes  the  final  method  in  which,  the 
rotor  geometry  specification  is  analyzed  for  control  of  the 
generated  voltage  wave  shape.  The  reluctance  of  the  air 
gap  is  specified  as  the  traveling  rotor  is  positioned  under 
consecutive  stator  slots.  .lust  as  in  Section  VIII,  reluc¬ 
tance  versus  time  may  also  he  related  to  air  gap  distance 
versus  time  and  thus  to  rotor  geometry.  The  specification 
j  progresses  with  time  along  with,  the  required  wave  shape. 

Therefore,  individual  harmonic  components  arc  not  of  concern. 

MMh  to  flux  density 

The  total  VM1:  versus  time  which  was  developed  in  Section 
VJ I  is  a  useful  starting  place  for  this  development.  Tie 
reluctance  must  he  specified  so  that  the  resulting  magnetic 
flux  density  will  he  in  accordance  with  the  magnetic  flux 
density  required,  as  was  discussed  In  Section  III.  Air  gap 
dominance,  the  lack  of  fringing,  and  the  finite  permeance 
effects  arc  assumed  here  as  they  were  in  Section  V ! !  I  . 

The  reluctance  at  any  slots  is  the  same  function  as  the 
reluctance  at  the  previous  slot  except  for  a  time  delay  of  At. 
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Therefore,  the  total  flux  density  may  be  expressed  as: 


where 


BT0T(t) 
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n  is  the  individual  stator  slot  number 
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Reluctance  Specification 

When  the  gap  between  salient  poles  is  under  a  slot,  the 
reluctance  at  that  slot  may  be  approximated  as  infinite.  When 
the  rotor  turns  so  that  its  edge  is  between  the  first  and 
second  conductor  of  the  relevant  coil  group.  Equation  (56) 
has  a  non-zero  term  only  for  n=l.  The  desired  increment  of 
flux  density,  AB,  and  the  available  MMF  determine  the  reluc¬ 
tance  required  while  the  rotor  edge  is  between  slots  1  and  2. 

MMF  (t) 

Rl(t)  =  — bTFT  (57) 

After  a  delay  of  t,  that  same  reluctance  will  appear 
under  the  second  slot.  The  MMF^  at  (t+At)  and  the  already 
specified  reluctance  will  determine  the  flux  density  at  the 
second  slot.  The  remainder  of  the  desired  total  flux  density 
must  come  under  the  first  slot.  By  using  MMF(t+At)  and  the 
residual  requirement  of  reluctance,  AB,  in  Equation  (57),  the 
reluctance  at  the  next  increment  may  be  determined.  The 
reluctance  specification  process  is  repeated  as  the  rotor 
progesses  under  the  coil  group  and  is  completed  once  the  new 
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rotor  position  is  tt  radians  from  whore  it  was  when  the 
specification  began. 

Ihe  specification  of  the  reluctance  using  this  method  is 
very  simple.  However,  insight  is  required  in  the  choice  of 
the  phase  voltage  wave  shape.  After  the  trailing  edge  of  the 
rotor  pole  face  passes  the  first  slot,  the  reluctance  specifi¬ 
cation  is  complete.  The  phase  voltage  and  the  required  mag¬ 
netic  flux  density  must  come  from  the  remaining  MM!  time 
variation  and  the  reluctance  specified  to  this  point. 

tixamp  1  e 

See  Appendix  I  f Generator  Design)  for  a  reference  on  the 
generator  spec i f icat ions  used  in  this  example. 

The  proof  that  not  just  any  phase  voltage  wave  shape  will 
work  is  seen  by  examination  of  the  example  developed  so  far. 

Table  5  gives  the  time  variation  of  the  total  MMF  at  each  slot 
for  the  Figure  5  voltage  variation.  At  w/t  =  -5r/15,  none  of 
the  pole  arc  is  under  the  phase  C  coil  group.  At  et =  -  4a / 1 5  , 
the  pole  will  be  under  the  first  conductor  of  the  coil  group. 

The  MMF  at  wt  —  Tir/lS  is  negative.  Additional  rotor  MMF 
could  be  added  here  in  order  to  ensure  a  positive  magnetic 
flux  density  with  a  realizable  reluctance.  However,  there  is 
another  difficulty.  The  desired  magnetic  flux  density  and 
therefore  the  desired  voltage  wave  shape  is  flat  between  .  t  --  -  2  /  1  5 
and  wt  =  +  .v;  /  1  5  .  At  tot  =  -2ir/15,  a  certain  combination  of  reluc¬ 
tances  will  be  under  slots  1,  2,  and  A  and  with  t he  MMFs 
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(Table  5  entries)  will  produce  the  desired  magnetic  flux 
density.  However,  at  wt=-Tr/15,  the  reluctances  are  exposed 
to  greater  MMFs.  The  new  increment  under  slot  1  also  sees 
a  positive  MMF.  Then  the  reluctances  which  provide  the  correct 
amount  of  flux  density  at  out=-2Tr/lS  are  bound  to  provide  too 
much  at  cot=-Tr/15.  Therefore,  this  choice  of  a  phase  voltage 
wave  shape  is  not  consistent  for  wave  shape  control  using 
only  rotor  geometry. 

Triangular  Phase  Voltage  Wave  Shape  Example 

The  flexibility  of  the  phase  voltage  wave  shape  for  a 
required  load  voltage  wave  shape  as  discussed  in  Section  III 
is  valuable.  Figure  9  shows  a  phase  voltage  wave  shape  which 
produces  a  DC  load  voltage.  Notice  that  it  has  an  upward 
ramp  i.ng  voltage  as  the  number  of  slots  with  finite  reluctances 
increases  and  a  downward  ramping  voltage  as  the  number  of  slots 
with  finite  reluctances  decreases.  In  this  example,  the  peak 
phase  voltage  equals  the  load  voltage  and  the  pole  arc  is  only 
tt / 3  electrical  radians.  Tables  7,  8,  and  9  display  the  voltage, 
the  magnetic  flux  density,  the  current,  and  the  MMF  variation 
at  each  slot  versus  time. 

The  reluctance  for  this  new  wave  shape  is  introduced 
under  the  first  conductor  at  <jjt=  -  3tt  /  30  .  The  rotor  pole  is 
centered  under  phase  C  at  a>t  =  -5r/30.  It  is  under  only  one 
conductor  at  u>t-13ir/30  as  it  is  leaving  the  coil  group. 


61 


Table  10  shows  the  progressive  calculation  of  the  reluc¬ 
tances  required  for  the  Figure  9a  triangular  voltage  wave  shape. 
Notice  that  Table  10  stops  with  the  determination  that  a  nega¬ 
tive  reluctance  is  required.  A  rotor  MMF  of  1400  ampere-turns 
was  chosen  instead  of  the  533  ampere-turns  used  previously  in 
order  to  prevent  the  calculation  of  a  negative  reluctance  for 
the  very  first  entry.  This  wave  shape  choice  is  not  appropriate 
for  rotor  geometry  control. 

From  the  shaded  area  of  Table  9,  it  is  observed  that  an 
equation  for  the  magnetic  flux  density  (similar  to  F:quation 
39)  may  be  written  for  each  line  (at  each  time  increment). 

Nine  equations  with  the  specified  MMF  determine  the  required 
reluctances.  Yet,  only  five  reluctances  are  variable.  The 
armature  current  and  rotor  MMF  provide  two  additional  variables. 
It  is  clear  that  design  using  the  method  of  this  section  is  a 
form  of  over- spec i f icat ion .  There  might  be  a  waveform  that 
fits  all  of  the  equations  and  provides  the  desired  load  vol¬ 
tage,  but  none  were  found. 


Table  7  Normalized  Phase  Current  and  Voltage 
(Triangle  Wave  Shape)  (Figure  9a) 
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(wt) 
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B 

_ 
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Table  8  Normalized  Armature  Reaction  MMF  Phase 
C  Slots  (Triangle  U'avre  Shape) 
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Table  9  Total  MMF  at  Phase  C  Slots  (Triangle  Wave  Shape) 
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Figure  9  Triangle  l’liase  Voltage  and  Current  Wave 
S  h  a  p  e  s 


X .  Other  Magnetic  Effects 


This  section  discusses  the  effects  of  flux  leakage, 
flux  fringing,  and  flux  saturation  on  the  desired  generated 
voltage  wave  shape.  Finite  Element  analysis,  an  important 
analytical  tool,  is  introduced  as  a  computer-aided  design 
method  used  in  fine  tuning  and  checking  the  generator 
design.  One  example  of  a  magnetic  finite  element  analysis 
program  is  introduced. 

Introduction 

All  of  the  procedures,  methods,  and  examples  in  this 
thesis  effort  assume  zero  flux  leakage,  no  flux  fringing, 
and  no  saturation  in  the  generator's  magnetic  structure. 
These  are  important  assumptions  because  calculating  flux 
leakage,  flux  fringing,  and  saturation  require  a  knowledge 
of  the  complete  magnetic  structure;  something  that  is  not 
available  in  the  initial  design  stages.  Generator  design 
is  by  nature  an  iterative  process  due  to  these  three 
problems . 

When  a  generator  is  designed  for  a  specific  load,  the 
system  inductance  is  an  important  specification.  The  svs 
tern  inductance  depends  on  the  load  inductance,  the  trans¬ 
mission  line  inductance,  and  also  the  generator  inductance. 
The  generator  inductance  can  easily  be  the  most  significant 
of  the  three  inductances.  The  generator  inductance  depends 


on  the  geometry  of  the  magnetic  circuit.  Designing  the 
magnetic  circuit,  however,  depends  on  a  knowledge  of  the 


system  inductance.  Normally,  a  "first  cut"  design  is  made 
from  which  a  "first  cut"  system  inductance  can  he  calcu¬ 
lated.  Using  this  inductance,  a  "second  cut"  design  is 
made  from  which  a  "second  cut"  system  inductance  can  he 
calculated.  This  process  continues  until  it  converges  on 
a  final  design. 

This  thesis  effort  is  concerned  only  with  the  "first 
cut"  design.  The  rest  of  the  design  procedure  is  left  for 
further  study. 

Leakage  Flux 

In  the  generator,  the  rotor  magnetic  flux  is  generated 
in  the  rotor  coils.  The  flux  travels  through  the  rotor 
poles  and  into  the  stator  magnetic  structure  cutting  the 
armature  coils  and  thus,  generating  the  output  voltage. 

In  a  salient  pole  generator,  there  is  a  large  air  gap 
between  two  adjacent  poles.  The  magnetic  reluctance  is 
very  high,  but  '  ^s  have  a  finite  value.  Since  there  is 
an  MMF  drop  betwc  two  adjacent  poles,  there  will  be  mag¬ 
netic  flux  crossing  the  gap  between  the  poles  bypassing 
the  stator.  All  of  this  leakage  flux  which  bypasses  the 
stator  does  not  have  a  significant  impact  on  the  shape  of 
the  desired  output  voltage  waveform,  although  it  does 
increase  the  geneiator  inductance,  X  .  To  compensate  for 


this  leakage  flux,  all  the  designer  needs  to  do  is  increase 
the  number  of  ampere  turns  in  the  rotor  coils. 

Fringing  Flux 

Another  important  assumption  made  in  this  thesis 
effort  is  that  all  the  flux  generated  in  the  rotor  passes 
through  the  air  gap,  past  the  armature  coils  and  into  the 
stator  core  in  a  totally  radial  direction.  This  is  an 
important  "first  cut"  assumption  and  can  be  corrected  for 
in  further  design  iterations. 

In  a  realistic  generator,  the  majority  of  the  magnetic 
flux  does  go  in  the  radial  direction.  However,  there  is 
a  significant  amount  of  flux  which  does  not.  This  flux 
crosses  the  air  gap  in  directions  defined  by  the  magnetic 
flux  density  vector  ?  causing  a  distortion  of  the  desired 
output  voltage  wave  shape.  The  flux  can  also  short  cut 
across  several  armature  slots  instead  of  going  straight 
through . 

Accounting  and  correcting  for  this  fringing  flux  so 
that  the  correct  flux  links  the  armature  coils  is  no  simple 
matter.  Where  and  how  the  flux  fringes  depends  on  the 
shape  of  the  generator's  magnetic  structure.  The  correc¬ 
tion  for  the  flux  fringing  is  therefore  an  iterative  design 


process . 


Magnetic  Saturation 

The  last  important  assumption  made  in  this  thesis 
effort  is  that  the  relationship  between  the  MMF,  the 
reluctance  of  the  magnetic  material,  and  the  resultant 
flux  density  in  the  magnetic  material  is  linear.  This 
last  assumption  said  another  way  means  that  the  analysis 
in  this  thesis  effort  ignores  the  effects  of  saturation  of 
the  magnetic  material. 

In  a  realistic  generator,  this  relationship  is  defi¬ 
nitely  non-linear.  When  the  MMF  dropped  across  the  magne¬ 
tic  structure  is  a  given  value,  a  certain  amount  of 
magnetic  flux  will  "flow'"  through  the  material.  As  the 
MMF  increases,  the  flux  flowing  through  the  material  also 
tends  to  increase.  However,  at  high  enough  values  of  flux 
density,  the  magnetic  material  saturates  and  more  and  more 
MMF  is  required  for  smaller  increases  in  the  flux. 

As  was  true  for  the  flux  fringing,  the  locations  and 
the  severity  of  the  flux  saturation  are  determined  by  the 
structure  of  the  magnetic  circuit.  Compensating  for  this 
problem  can  also  be  an  iterative  design  procedure. 

Finite  Element  Analysis 

The  calculations  involved  in  designing  a  generator 
can  be  exceedingly  complex  and  cumbersome.  This  fact  makes 
the  digital  computer  a  welcome  design  aid.  One  particular 
procedure  ideally  suited  for  computer-aided  design  is 
Finite  Clement  analysis. 
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In  Finite  Element  analysis,  a  cross  section  of  the 
generator  is  taken  perpendicular  to  the  axle.  Then  the 
assumption  is  made  that  any  leakage  flux  or  fringing  flux 
is  in  the  plane  of  the  cross  section.  That  is,  there  is 
no  leakage  flux  or  fringing  flux  between  any  two  "adjacent, 
infinitely  thin"  cross  sections.  The  two  dimensional  cross 
section  is  then  divided  up  into  triangular  finite  elements. 
Each  element  is  restricted  to  one  particular  type  of  mate¬ 
rial  (steel,  air,  insulation,  ...).  During  the  computa¬ 
tions,  each  triangular  finite  element  is  considered  to 
have  a  constant  magnetic  flux  density  over  the  entire 
element.  See  (Ref  8)  and  (Ref  13)  for  further  information 
on  the  Finite  Element  method. 

MAGNETIC  is  a  computer  program  that  uses  the  Finite 
Element  analysis  technique  (Ref  4)  to  analyze  magnetic 
structures.  MAGNETIC  computes  the  magnetic  fields,  induced 
currents  and  losses,  and  the  forces  of  a  magnetic  device 
that  has  been  represented  as  a  collection  of  finite  ele¬ 
ments  . 


"MAGNETIC  aids  the  design  of  magnetic 
devices  because  accurate  prediction  of  devices 
performance  usually  requires  knowledge  of  the 
magnetic  field.  MAGNETIC  is  a  means  of  calcu¬ 
lating  the  field  that  allows  for  fringing, 
saturation,  permanent  magnets,  and  complicated 
shapes.  Because  the  program  includes  non-linear 
B-H  curves  and  computes  inductances,  forces,  and 
eddy  current  losses,  it  can  help  the  designer  to 
minimize  the  materials  used  in  a  device  and  maxi¬ 
mize  its  efficiency."  (Ref  4:11) 


Two  problems  in  generator  design,  flux  fringing  and 
saturation,  force  the  design  procedure  to  be  iterative. 
MAGNETIC  can  be  a  welcome  aid  in  correcting  for  flux 
fringing  and  saturation  as  well  as  flux  leakage.  After 
the  designer  has  completed  the  "first  cut"  design,  the 
design  can  be  analyzed  by  MAGNETIC.  MAGNETIC  shows  through 
its  output,  where  and  to  what  extent  flux  fringing,  satura¬ 
tion,  and  flux  leakage  are  occurring.  With  this  informa¬ 
tion,  the  designer  is  able  to  take  corrective  action  and 
come  up  with  a  "second  cut"  design.  This  procedure  con¬ 
tinues  until  an  acceptable  design  is  obtained. 

In  addition  to  tabular  listings  of  the  various  output 
parameters,  MAGNETIC  produces  a  Mesh  plot  and  a  Flux  plot. 
The  Mesh  plot  portrays  the  selected  generator  cross  section 
showing  all  the  triangular  finite  elements.  The  Flux  plot 
shows  the  lines  of  magnetic  flux  in  the  generator  cross 
section.  Figure  10  is  the  Mesh  plot  for  the  sample  gener¬ 


ator  [Reference  Appendix  I).  The  Flux  plot  of  the  sample 
generator  is  not  shown  because  that  portion  of  MAGNETIC 
currently  does  not  work. 


XI.  Conclusions  and  Recommendations 


Conclus ions 

The  purpose  of  this  thesis  effort  was  to  investigate 
the  use  of  rotor  geometry  as  a  specification  in  the  control 
of  generator  voltage  and  flux  density  harmonics.  It  was 
found  that  armature  reaction  MMF  could  be  expressed  as  an 
infinite  set  of  traveling  waves  with  an  infinite  set  of 
velocities  (Section  IV)  and  (Hquation  23).  Because  the 
rotor  MMF  has  only  the  velocity  of  the  rotor,  it  cannot 
cancel  armature  reaction  everywhere  in  the  generator 
throughout  all  time  (Section  VII). 

Rotor  geometry  was  put  to  use  in  a  limited  control  of 
flux  density  harmonics  at  the  stator  slots.  The  reflection 
of  the  MMFs  present  at  all  of  the  slots  of  a  coil  group 
into  a  single  slot  allows  the  possible  calculation  of  an 
appropriate  rotor  geometry.  This  reflection  turned  out  to 
exceed  the  resources  of  this  thesis  effort  and  was  not 
accomplished  (Section  VIII). 

Specifying  the  air  gap  reluctance  incrementally  as 
the  rotor  moves  under  successive  stator  slots  was  also 
tried  as  a  means  of  harmonic  control.  This  method  was 
unworkable  for  all  wave  shapes  chosen  because  it  produced 
negative  MMF  requirements  and  conflicting  geometry  spcci- 
f icat  ions . 
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Recommendations 


The  control  of  flux  density  is  an  important  consid¬ 
eration  when  working  to  improve  generator  efficiency  and 
performance.  Because  of  the  potential  that  flux  density 
control  has,  several  recommendations  for  further  study 
are  presented. 

Before  any  hardware  is  fabricated,  it  is  recommended 
that  the  designer  calculate  the  benefits  due  to  the  flux 
density  control  design.  Knowing  that  energy  is  conserved 
between  the  input  and  the  output  of  the  generator,  the 
designer  will  be  able  to  calculate  the  design's  value. 

A  recommended  approach  for  reduction  of  losses  is  to 
investigate  flux  density  control  at  the  damper  bars  instead 
of  the  stator  windings.  This  investigation  however,  will 
not  control  the  generator  output  harmonics. 

A  different  approach  might  be  to  use  stator  geometry 
as  a  design  specification  in  harmonic  control.  This  will 
not  control  the  MMF  at  all  time  and  all  space  by  itself, 
but  may  provide  easier  control  of  the  harmonics  at  a  given 
location . 

The  method  in  Section  VII  could  be  further  investi¬ 
gated  if  a  method  were  derived  for  the  division  of  one 
infinite  series  with  progressive  harmonic  arguments  by 
another  such  series.  The  Section  VIII  approach  could  be 
further  investigated  in  terms  of  a  general  wave  shape 


(see  the  heading  "The  Contradiction"  in  Section  VIII). 
hither  method  might  provide  a  workable  wave  shape  or  prove 
that  the  rotor  geometry  concept  is  in  general  impossible. 

The  general  requirement  for  phase  voltage  wave  shapes 
that  produce  a  l>C  output  from  a  three  phase  full  wave 
bridge  rectifier  is: 

T 

Voltage  (t)  +  Voltage  {-y-  +  t 1  =  a  constant  (58 
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s  in  tie 

a i r  gap  . 

The  final  recommendation  is  to  include  the  flux 
fringing,  flux  leakage,  and  magnetic  saturation  problems 
in  the  calculation  (only  if  the  calculations  are  successful 
without  these  three  problems)  for  the  desired  phase  voltage 
wave  shape.  It  should  be  recognized  that  the  problem  will 
now  be  iterative  and  Finite  Element  analysis  (Section  X) 
will  be  helpful . 
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Appendix  I  Sample  Generator 


Specifications  and  Dimensions  for 
the  Sample  Generator 
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The  development  of  the  rotor  geometry  modification 
example  in  this  thesis  requires  a  constant  air  gap,  salient 
pole,  synchronous  generator  to  begin  the  investigation. 

The  generator  selected  is  basically  the  one  designed  for 
the  final  exam  in  Air  Force  Institute  of  Technology  course 
EF  6.07  (Electromagnetic  Components  Design,  Spring,  1980). 

A  few  modifications  are  made  due  to  thesis  requirements. 

The  design  procedure  for  the  sample  generator  is  not  pre¬ 
sented.  Specifications  for  the  sample  generator  are 
presented  in  Table  11.  A  cross  sectional  drawing  (per¬ 
pendicular  to  the  axle)  of  the  generator  is  presented  in 
Figure  11. 

Th "  design  modifications  developed  in  this  thesis 
apply  equally  as  well  to  any  salient  pole,  synchronous 
generator  as  they  do  in  the  generator  used  to  begin  the 
invest igat ion . 

Notice  that  the  generator  design  is  in  inches.  \11 
units  shown  in  this  appendix  are  converted  to  the  MKS 
system  for  use  in  the  body  of  this  thesis. 
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Table  11 


Sample  Generator  Dimensions  and  Specifications 


Machine  Specifications 

WYE  connected  three  phase  output 

460  volts  RMS  line  to  neutral  phase  voltage 

72.464  amperes  RMS  phase  current 

Load  power:  100  kilowatts 

Frequency:  400  hertz 

No  rotor  or  stator  vents 

Straight  slot  stator  design 

60  slots 

4  poles 

0.015  inch  constant  air  gap  over  the  poles 

Synchronous  speed:  12000  revolutions  per  minute 

Rotor  tip  speed:  393  feet  per  second 

Rotor  MMF :  533  ampere  turns  per  pole 

Maximum  flux  density:  1  Tesla 

Pole  arc:  60  mechanical  degrees,  a  =  0.667 

Rotor  length:  25.347  inches  P 

Trancor  "T"  steel  with  0.007  inch  laminations 

Synchronous  salient  pole  design 


Winding  Specifications 

Six  internally  generated  phases,  m  =  6 
Armature  conductors  have  no  skew,  k'sv  =  1.0 
q  =  5 

One  series  current  path  per  output  phase 

Full  pitch  stator  winding,.  W/T  =  1.0,  K  =  +  1.0 

Two  layer  lap  design  ™  ™ 

0.107  by  0.107  inch  square  conductors 

15000  amperes  per  square  inch  maximum  current  density 

80  volts  per  mil  minimum  insulation  thickness 

T  =15  slots 
P 


Assumptions  (for  "first  cut"  calculations) 


Zero  flux  leakage 
Zero  leakage  reactance 
PF  =  1.0 
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Figure  11  Cross  Section  of  the  Sample  Generator 
(Specifications  and  Dimensions) 
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Appendix  I 


RGRID  and  TRIA  Card  Generation 


Computer  Program  to  Generate  RGRID  and  TRIA  Cards 
Required  by  MAGNETIC  (Ref  4)  for  the  Air  Gap  and  Pol 
Shoe  Regions  of  the  Sample  Generator  (Appendix  I) 


Appendix  II  RGRIP  and  TRIA  Card  Generation 


The  following  program  listing  generates  the  RGRID  and 
TRIA  cards  required  by  MAGNETIC  (Ref  4)  for  the  air  gap 
and  pole  shoe  regions  of  the  sample  generator  (Appendix  I). 

Refer  to  Figure  12  and  to  the  program  listing  for  an 
explanation  of  how  to  use  the  program.  Basically,  to 
modify  the  program  for  another  generator,  the  following 
must  be  accomplished: 

1)  Fill  in  the  appropriate  numbers  in  the  three  data 
statements . 

2)  Calculate  values  for  the  theta  and  radial  coordi¬ 
nates  for  the  points  A,  B,  C,  and  D. 

3)  Write  a  statement  function  that  describes  the 
surface  of  the  pole  shoe. 

The  internal  program  structure  is  not  discussed  because 
it  is  not  important  for  the  program's  use. 


86 


o£  /—v 

uuo 

CO  JW 
S  <  H 
3  U  < 
Z  CO  cc 
CU 

Q  O  O 
mHU 
a £  < 

Of-X 
cu  o  cu 
z 

<  > 

CO  .J 
t/5  m  H 
W  < 

hOW 
<  Z  cu 

CJ  >-H  u 
►-I  3c 
Q  <  Cl 

z  cu  < 

I-H  Q  O 


1  CO  CU 

r-4 

.  •  t— (  1 — 1 

/ - \ 

co 

CO  CTlIZ  < 

< 

cu 

CU  <ni|h  w 

E— 

H 

CU 

t- 

O  t — >, — > 

T 

z 

Zhn 

H 

cu 

' — ^ 

s: 

u 

CL 

cu 

OS 

co 

3 

CO 

\ 

CL,  CO 

c«s  o 

H 

E~  < 

Z  H 

CU  CO 

S  I 

u 

cu  I 

co  I 


CL 

< 

C J 

05 

QC 

CL 

CU 

< 

H 

CU 

Z 

z 

cu 

CU 

H 

u 

s 

< 

o 

CL 

CL 

cu 

O 

CU 

co 

CU 

v 

►u 

o 

Q 

Cl 

O 

i 

\  cn 
E— •  OQ  O-J  I 
Z  CL 
CU  LU 


CU 

CU  O 


Figure  12  Cross  Section  of  the  Sample  Generator  (RGRID  and  TRIA  Cards) 


o  ii  (.1  n  ('  n  o  ft  tt  (i  t'  tt  it  it  it  <t  ct  o  it  o  o  r>  t'  O  C'  ti  o  o  o  «i  O  n  o  n  o  o  it  it  it  ct  ooro  o  tt 


PROGRAM  THESIS1  ( INPUT, OUTPUT , PUNCH,  T  A  PEE  =  OUTPUT  »T  APP?=PLHCH) 

TUTS  PROGRAM  GENERATES  ONLY  THE  A IR  SA*  AND  30LE  SHOE  RGRTO 
AK-0  CAPOS  REQUIRED  BY  THE  PROGRAM  MAGNETIC  (REF  to  FOR 

juc  SAM°L  S  GENERATOR 

VARIABLE  DImENSI OSS 

ang  -  angle  for  individual  rgriq  cards 

XEY  -  KEYri  DESIGNATES  A  TRIA  IS  AIR,  KEY  =  "  DESIGNATES  rT EEL 
L  MG  -  LNG  +  i  =  FIRST  GRID  NUMBER  ASSIGNED 

INS  -  LNS+1  =  LOWEST  G;*ID  NUM^Er.  ON  THE  STATOP  AIR  G*P  BOUNDARY 

LMT  -  LNT+l  =  FI " ST  TkIA  NUMBER  ASSIGNED 

IPG  -  LPGM  =  FIE  ST  GRIO  NUMBER  ON  THE  RB  ARC 

nf  -  MUMPER  OF  c-GRIO  °C INTS  ON  SEGMENT  KF  (INCLUDES  THE  ENDS) 

Nr A  -  NUMBER  OF  PGRID  POINTS  ON  SEGMENT  KFA  ( INCLUDES  THE  ENOS) 

Mtrg  _  NUMBEV  OF  RGRID  POINTS  ON  SEGMENT  KFB  (INCLUDES  THE  ENDS) 

NG  -  RGRIO  NUMBER  TO  BE  PUNCHED 

NGi  -  NUMBER  OF  PGRID  POINTS  ON  SEGMENT  KOI  (INCLUDES  THE  ENDS) 

T‘P  -  N1.JMBPF  CF  FORTE  POINTS  ON  SEGMENT  KP  (INCLUDES  THF  ENDS) 

MPA  -  NIJMBE-  OF  RGRID  POINTS  ON  SEGMENT  KP&  (INCLUDES  THE  ENDS) 

NP3  -  NUMDEF  OF  RGRID  POINTS  ON  SEGMENT  KP*  (INCLUDES  THE  ENDS) 

NT  tp  I A  HUMBER  TO  BE  PUNCHEO 

rI  -  CONSTANT  3 . 141593653 5  ... 

RAO  -  RADIAL  DISTANCE  TO  THE  POLE  FACE  SURFACE 

BG1  -  RADIAL  DISTANCE  TO  THE  MI DDL  E  OF  THE  AIR  GAP 

YAB,  XCC,  XCD,  XT,  XY,  X7,  AND  YCC  ARE  USED  TO  CALCULATE 
THE  COORDINATES  FOR  THE  POINTS  A,  9,  0,  AND  D. 

PA,  PB,  RC,  AND  PC  ARE  THE  RADI  •"  l  COORDINATES  OF  A,  p,  C,  AND  D. 
TA,  TB,  TC,  AND  TD  ARE  THE  THETA  COORDINATES  OF  A,  B,  C,  AND  Q. 

THE  FOLLOWING  VARIABLES  ARE  USED  INTERNALLY  TN  THE  PROGRAM  ONLY  t 
DKTB  DNFP  ONFPA  ONFPB  DNGFM?  DNSGi 

I  KB  KT  NFAM1  NFBHi  NFM1 

NG1M1  NPAM1  NPBMi  N»M1 

THE  FOLLOWING  ARE  DESIGNATED  LINE  SEGMENTS  -  SEE  FIGURE  12 

KF  KFA  KFB  K31  KP  <PA 

K  p  B 

SUBROUTINES  F EQUIP  ED J 

TRIA1  -  PUNCHES  TRIA  CARDS  FOP  rv»E  1  FINITE  ELEMENTS 
TRIA  2  -  FUNCHES  TRIA  CARDS  POR  TYPE  2  FINITE  ElEMFNTS 
SEE  T RI A  1  AND  T K I A 2  FOR  FURTHER  EXPLANATION. 

DOUBLE  PRECISION  TA, T D, TC  , 7D , RA , R3 , RC , PD 

OOU3LE  PRECISION  ANG , RA  D , PI 

DOUBLE  PRECISION  XY, X Z, XT  , XCC , YCC, XCD , XA 3 
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O'TA  NFA,NF,NF°,NPA,NP,NP9,NG1,NS  1 1 1 ,  *  1 , 1 1 ,  f  ,  3 3  ,6  ,61, 61/ 
TATA  LNG , LNT  ,  LCG ,  LNS  /275,4C >,23 , 214/ 

DATA  PGl  /3 


0  juts  STAT  FMc  KT  PUNCT  ION  REPRESENTS  THE  SURFACE  OF  t  Mr  POLE 
:  S^E.  IT  IS  iHF  DISTANCE  VERSUS  X  (HE  ORIGIN  IS  IN  T Hr 
C  C-NTH  OF  THE  POLE)  A  90V  E  THE  PB  ARC  ♦  IN  'T  HIS  CAST,  t  CONSTANT 
C  AT»  G«o  has  PEEK  SPECIFIED.  THE  ONLY  RESTRICTION  13  THA T  THE 
C  °CLE  PACE  hay  NOT  HAVE  ANY  PORTIONS  HAT  A  RE  CONSTANT  TH  FT  A 
C  LTNES. 

SURF < THETA )=• . 21 2257 075 4  7377 1*- 1 7 E <*7 5236Dr 
t«''  FORMAT  (•‘I**) 

IV1  FORMA  T  (”  M,  "KGr  ID”  »5X»I3»2(5X»C15«1D  ) 

1  ’?■*»  FORM AT  ("r. Gkl 0", 5 X  ,  13 , 2(  5X  ,F1  O.iv  )  ) 

GTVEN  THE  EQUATION  FOR  THE  POLE  SHOE  SURFACE,  CALCULATE  THE 
COORDINATES  FO^  THF  POINTS  A,  9,  C,  AND  D. 

THE  CONSTRAINTS  A R E I 

1)  THE  ARC  LENGTH  BETWEEN  POINTS  A  AND  2  ^UST  BE  6C 
MECHANICAL  DEGREES. 

2)  THE  LINE  Sr  GHENT S  A9  AND  DC  MJST  PE  PARALLEL  TO  A 
LINE  THROUGH  THE  CENTER  OF  THE  GENERATOR  DIVIDING 
THE  POLE  IN  HALE, 

PI  =3.  1415926 63E  8979328346264335327950  0 

TA  =  7E.•’D<-■‘PI/leJ.'!D:, 

T3  =  lc.  ’0  'PI/IE?. IDO 
XY  =  OSORT(  j. tO.  *  3. 3601*3. 3600 
X?=OSORT (j.EOf*  1.3 200*1 #  32DC ) 

XT  =  QS  ORT  (  c  .  *  0>  *  0.1 95  03*6.19  EDO) 

YCC=XY+X7+XT 
YCC=X7-XT 

RC-DSQRT ( ( XCC* XCC) ♦  <  YCC*YCC>  ) 

RB  =  RC 

PA=SURF(TA>  +RB 
RO=SURF (TO) +RP 

XA9=RA*OSIN(3i.aOl*PI/18u.  '0  C) 
jO=(A5.CDw*FI/13.}.CDi.)+DASIN  (XA9/RB) 

XC 0=P  0*0SIN  (3;-  ,  I OO’PI /1 8 0  •  co  ;> 

TCs(4  5,l.DC*PI/13r  .CDu)-OASIN(XCVRC) 

WRITE  (6,1  VIC) 

PUNCH  THE  RGRIP  CARDS. 

RGRID  POINTS  APE  TO  PE  EQUALLY  SDACED  ON  GIVEN  SEGMENTS. 

REMOVE  IMF  C  IN  COLUMN  1  PRIOR  10  ALL  WRI TE( 7 ,  ...  STATEMENTS 
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3  If  PUNCHED  CARO  OUTPUT  IS  DESIRED. 

r* 

Z  PUNCH  R3CTD  CAFDS  FOR  SEGMENT  K3  R  -  ONLY  THE  ?  I G  HT  H  A  N  D  ENT  POINT 

/* 

U 

ng=l°g 

DO  1;  ]=i,NPBMi 

ANG=  (  TO/*  P- '■•1)  ’  <I-l>Mi8j.  D  /PI) 

NG=NG+i 

W  RIT  c (m,1  1 r  )  VG ,RB, ANG 
0  WRITE  (7,1  u2(  )  KG,R8,ANG 

tr  CONTINUE 

Z  o(!NCH  ?G  R 1 0  CAf-CS  FOR  SEGMENT  <3  -  ON.Y  THE  RIGHTHANO  END  =>Oi:JT 


N°H1=  N°-l 

DO  2'.  I  =  l,NPMi 

ftNG=  (TC*(  (  (TB-7C)  /NPMt)  *  (1-1)  ))  M  i30.uD’j  /  PI) 

NG=NG  *1 

WRITE  (S,l.'lr)  NG,P9,  ANG 
0  WRIT  c  (7*  1«'  2f  )  NG,RE,ANG 

?'  CONTINUE 

C  PUNCH  ,RG p 1 0  CARDS  FOR  SEGMENT  KPA  -  90TW  END  POINTS 

•J 

N DA  HI  r  NPA  -1 

DO  3C  1  =  1, NOA 

ANG=(T3+(  ((  (PI/2.  !!00)-T8)/NPAM1)  *  (1-1 ) ) )  »  (  18  C.!  DP/PI) 

NG=Nr-  *-1 

WRITE  (5,1 ilO)  NG,R5,ANG 
3  WRITE (7, 1*2' )  NG , RB, ANG 

71  CONTINUE 

•* 

Z  PUNCH  RGRID  CARDS  FOR  SEGMENT  ^9  -  ONLY  THE  R1GHTHAND  END  DOI NT 
C 

NG=L  NG 
NRRH1=NFB-1 

DO  Ar  1=1, NPBM1 

ANG= (TD/NFBM1) r  ( I - 1) ♦ ( 1 8 0 .  CO . /°I) 

NG=NG+1 

WRITE  (E>,  lLlf )  NG  ,  RD,  ANG 
C  WRITE  (7,li-2f  )  NG  ,  RD,  ANG 
k ?  CONTINUE 

r» 

C  PM NCH  R3RI0  CARCS  for  SEGMENT  <-  -  DOTH  END  3DINTS 
NCM1=NF-1 

DO  51  I=1,NFM1 

ANG=  (TD-M  ((TA-TO)/NtrMl),»(I-l))),,(13J.CO:/:>n 
RAO=SUPF(ANGH-cB 
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nn  n  n  o  c*  <•>  r>  o  o  n  o  n  o  o  o  o  O  n  o  o  f> 


NG  =  NG  «-i 

WRITF(*S,1j1!:>  N'G  ,  RAO  t  ANG 
WRITE(7,1CZI>  N'G  »  RAD  » ANG 
~r  CONTINUE 


omnch  rs-id  ca»-.rs  for  segment  kpa  -  qdm  e^o  points 


Nr  AMi -NEA -1 

DO  6'  I=1,NFA 

ANG=  (T  !U<  (  (  (PI/2.  .CMJ-T  A)/MrAHi)  *  (  I  - 1 )  )  )  »  (  IP  .  Or  /PT) 

NG=MO  4-i 

WRIT"  (  S , 1 w 1 ’ )  NG,PA,ANG 
WRITE  C*,  lw2.  )  N'G,  RA,  ANG 
"  CONTINUE 


P"NCH  s>Gr-I0  CARDS  FOR  SEGMENT  <31  -  30TH  END  FOINTS 


NG1M1=NG1 -1 

DO  7.  I  =  1 » NG 1 

ANG=  (  (  (Pl/2.  .  D  )  /NG1MD*  (1-1 )  )•  (IS!  .  ,Or  /=>T  ) 

NG  =  NG  >1 

WRITE  ( o , 1 0 1 T  )  NG,RG1,ANG 
WRITE  (7,1.2'*)  EG  » RGi »  ANG 
71  CONTINUE 

ALL  Oc  THE  RGPIO  CARDS  APE  NOW  PUNCHED 
PROCEED  TO  PUNCH  THE  TRIA  CARDS 

RrwO  Mc  THE  C  IN  COLUMN  1  PRIOR  TO  ALL  WRIT  E.  ( T ,  ...  STATEMENTS 
IF  PUNCH  Cfl  OUTPUT  IS  DESIRED. 

PUNCH  TRIA  CARDS  FOR  REGION  T1  (DETWEEN  THE  OUADRATURF  AXIS 
AND  THE  RIGHT  SIDE  OF  THE  POLE  SHOE 

KEY=1 
NT  =L  N'T 
KT  =1 
KD  =  1 

DNFPP=cLOAT (NFP) /FLOAT (NPB) 

8"  OKTD=  FLOAT ( KT ) /FL 0 AT ( Kn ) 

IF  (HKTB. GT .PNFOP)  GO  TD  9c 

IF  (KT.EO.NFB)  GO  TO  IPO 

CALL  TRIAi  (KT,K3,LNG,LPG,NT ,KET) 

GO  TO  8T 

«<?  CALL  TRIA2  (KT , K B , LNG ,L PG ,NT , KEY ) 

GO  TO  8"» 

PUNCH  TRIA  CARDS  FOR  REGION  T2.  (THE  POLE  SHOE) 
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i  '  KEY=‘ 

L‘JG=L  NG+NFP-1 
L°G=LPG+NPP-1 
KT  =  1 
<3  =  1 

Q  JPe= rLC AT (KF) /FLOAT  (NP) 
i  1  ‘'<TO  =  cLOAT  (<T)/cL0AT  ( KP  ) 

Ic  (HKTP.GT.ONP3)  GO  TO  12 

; c  (YT.EQ.LF)  GO  TO  13' 

CALL  TPIfi  (KT  ,Kq»LNGfLPG,NiT  ,  KEY) 

GO  TO  Hi 

1?'  CALL  T  F  i  A  2  (KT,Ke,LNG,LPG,MT,KEY) 

GO  T0  il* 

°UNCH  TRI 3  CAPOS  FCR  " EGI ON  T3  (3ETWEEN  THE  QUAD FAT  URE  AXIS 
AMO  THE  LEFT  SITE  CF  THE  POLE  SHOE) 


l’*"  KEY=  1 

LNG=L  NG  +  NF-1 
L3G=  L  PG*1-  P-  1 
<T  =  1 
<1  =  1 

ONFPAsFLOAT (NFA) /float (N PA) 

1'  •’  0<T3=FLOAT  (<T) /  FLOAT  ( <9 ) 

Ic  (OKTE.GT.DNFPA)  GO  TO  15 j 

I f  (KT.EQ.NFA)  GO  TO  If, 

CALL  TPltl  (<T,K9,LNG,LPG,NT,<EY) 

GO  TO  U[ 

If  CALL  TRIA2  (KT » K? »LNG »L PG »NT  >KEY) 

GO  TO  if 

PUNCH  TRI  A  CARDS  FOR  REGION  T4.  (LOWER  HALF  OF  THE  AIR  GA  °) 

If"  KEY= 1 

L  NG=  L  NG+NFA 
L°G=LNG-NFA-NF-NF9+2 
<T  =  1 
<0  =  1 

DNGFM  2  =  FL  OAT  (NGD  /FLOAT  ( NFA +NF*HF R-2) 
l""’  DKT3=F|_0AT  <<T)  /FLOAT  (  <9  ) 

IP  (DKTP.  GT.DNC-F^)  GO  TO  13.' 

IF  (<T  .EC.NGi)  GO  TO  193 

CALL  TRIA1  (KTfKO,LNG,LPG,M,  KEY ) 

GO  TO  17 1 

19J  CALL  T 01 A 2  ( KT , < ? ,LN G ,L PG ,NT , <EY) 

GO  TO  171 

PUNCH  TRIA  CAROS  FOR  REGION  T3.  (UPPER  HALF  OF  THF  AIR  GA °) 
19 n  <EY=1 
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r 


LMG=LMS 

L°G=LPG+(NFA+SF+NFB-2) 

KTT  =  1 
<^=1 

r'’nGi=cLPA’  (NS) /FLOAT (NGl) 

?  '  0<TB=cL06T <*T) /FLOAT (<<?) 

lc  (nKTB.GT«3KrGl)  GO  TO  21 • 

I"  (KT.t iC . »-*S >  GO  TO  22' 

CALu  T3IA1  <m,0,LNG,LPGtNTf  <ry> 
GO  TO  2i  r 

?'  •  CALL  T 31 A 2  (K1  ,<P  ,LNG,LPG,K'T  ,><2Y) 
GO  TO  ?t  • 

2?f  STOP 
E.NO 
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J 

Z  TwiS  SUBROUTINE  PUNCHES  TRI  A  CAROS  FDR  TYPE  1  FINITE  ELEMENTS. 
C 

Z  S'-  TH'STSi  FOr  AN  EXPLANATION  DF  AL.  VARIABLES. 

Z  “-'MOVE  TH-  C  iv  COLUMN  1  »-:IO!  rU  Al»  WRITER,  ...STATEMENTS 
Z  I r  P'JNCHcO  OUTPUT  IS  OESIRfO. 

Z  Tver  j_  PTNITt  fc  L  E  f-  r  N  T  S  ARE  TRIANGLES  WITH  O0WNWSPP  POINTING 
C  A® “ XES  • 

J 

SUBROUTINE  TRiai  (  KT  ,KP,LNG,LPG,NT  ,  KEY  ) 

N  2  =  <" +LNG 
N1=N 2+1 
N2=KB  +LPO 
KT=KT  *1 
NT=N~ 

IF  (KEY, EC. i)  GO  TO  2° 

WRIT=:  <6,1.  1'  )  NT  ,  N1,N2,N3 
C  WRITE  (7,1-11)  v T , N 1 , N 2 , N 3 

1M‘  FORMAT  (**  "  ,  “T  ^  I A  “  j  -+  (  ;  X,  I  2  )  ) 

I'll  FORMAT ("TKI A", k (EX, 13)) 

F.ETIJP  N 

J 

Z'  WRIT  ■;  (G ,  1  •-  2  :  )  N’T  ,N1,N2,N3 
C  WRITE  (7,11'21)  NT  ,M,N2,N  ? 

1*2)  FORMAT  (**  " ,  *’TF  I A  “  »R  ( &  X,  13)  ,5X,'*AIR**) 

1"21  FORMAT(*,TRIA’*,t(EX,I3)  ,  5  X  ,"A  IR”) 

PET'JPN 

ENO 
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T  UIS  SU^~  OUT  1  N£  PUNCHES  TRIA  CAPUS  FOR  TTPE  ?  FINITP  ELFHENTS 


0  E~C  THESIS1  F  0-  AN  EXPLANATION  OF  A  _  L  YAR1A«_ES. 

•> 

0  S’-'-'V/P  c  tk  CfLUP.N  1  PRIOR  TO  ALL  •  .  .f  T  ^  T  EH  rNTS 

C  ’r  P'.INCHFO  OUTPUT  IS  RETIRED. 


0  "vrjE  ?  FTNITE  cLcNrNTS  ARE  TRIANGLES  WITH  UPWArT  PfTMT  ING 
C  Ac-x^.c:, 

SUBROUTINE  T R 1 A 2  (KT,I<9,LNG,LPG,NT,  KEY) 

Ni  =  KT  +LNG 
N2=<p+LPG 
Nf=N2+l 
<P  =  <P  +1 
NT=NT  +1 

IP  ( KEY»C0*  i)  GO  TO  2:1 

WRITE (5,1. li )  MT,N1,N2,N3 
0  WRITE  (”,1.11)  NT,N1,n2,N3 

1  M  '  F9RHATC*  ’* »  "T  - 1 A  ••  t »,  ( E  X ,  I  3 ) ) 

1 n  1  i  FORHA  T  <"TRIA**jA  ( 5  X  *  I  3) ) 

RETURN 

r\ 

J 

2°  WRIT'  («>,1.2C  )  NT,N1,N2,N3 

0  WRITE(7,i. 21)  NT »  Nil N2j  N 3 

1  ”  2r  F0P1AT("  ’*  i  *’T  ^  i  A  "  i  R  ( L»X|  I  3)  |  E  X  |  "A  I  *  ”) 

FORHAT  ("TP1  A’V  (rX,I3)  |5X,"AIR") 

RETURN 

ENO 
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because  the  rotor  MMF  has  only  the  velocity  of  the  rotor. 

The  MM  Is  present  at  all  t  lie  stator  slots  of  an  armature  coil 
group  are  reflected  into  a  single  slot  in  an  attempt  to  control 
the  air  gap  flux  density  harmonics  tie  selectively  positioning 
increments  of  rotor  reluctance.  This  method  fails  because  a  gen¬ 
eral  solution  for  the  reflected  flux  densities  could  not  he 
o  b  t  a  i  n  e  d  . 

Specifying  the  air  gap  reluctance  incrementally  as  the  rotor 
moves  under  successive  stator  slots  is  also  tried  as  a  means  of 
harmonic  control.  This  method  proves  unworkable  due  to  negative 
MMF  requirements  and  conflicting  geometry  specifications. 

Finite  Flemcnt  analysis  is  introduced  as  a  computer  aid  to 
the  design  of  a  magnetic  structure. 

Recommendations  for  further  studv  are  included. 
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